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1      INTRODUCTION 

1.1    GENERAL 

Today,  many  computer  programs  in  structural  analysis  are 
readily  available  and  each  is  designed  for  a  specific 
use.  Just  to  name  a  few  are  SAP90,  ETABS,  SAFE,  M- 
STRUDL,  MICROFEAP,  MULTIFRAME,  RISA-2D,  etc. . .  and  many 
more.  These  programs  are  generally  used  for  structural 
analysis  and  design.  Most  of  these  programs  have 
graphics  capability  which  provides  visual  presentation  on 
the  behavior  and  response  of  the  structure. 

While  most  of  these  programs  are  developed  using 
sophisticated  computer  methods  in  elastic  structural 
analysis,  few  are  designed  to  perform  elastic  second- 
order  analysis.  Its  use  can  be  of  considerable 
importance  in  the  design  of  various  structural  systems, 
particularly  frame  type  structures . 

In  this  report,  the  development  and  implementation 
of  computerized  elastic  second-order  analysis  of  frames 
is  presented.  The  technique  is  illustrated  by  several 
examples . 


1 . 2   OBJECTIVES 

The  primary  objective  of  this  report  is  to  develop  and 
implement  into  a  computer  program  the  concepts  of  elastic 
second-order  analysis .  An  additional  objective  is  to 
extend  these  concepts  to  perform  computerized  elastic 
stability  analysis . 

In  the  current  LRFD  Specification  for  steel 
buildings  '  1  ]  (1st  Ed.,AISC  1986),  Chapter  C,  Section 
CI  states  that,  "  second  order  (  P -Delta  )  effects  shall 
be  considered  in  the  design  of  frames.  "  Chapter  H, 
Section  2a,  of  the  Manual,  which  deals  with  beam-column 
design,  specifies  three  (3)  ways  to  determine  the 
required  flexural  strength  Mu  of  a  member: 

a.  M  u  may  be  determined  using  second-order 
elastic  analysis  using  factored  loads 

b.  M  u  may  be  determined  using  plastic 
analysis  that  satisfies  the  requirements  of  Sections  CI 
and  C2  of  Chapter  C  of  the  Manual .  Section  CI  deals  with 
the  second-order  analysis  while  Section  C2  deals  with 
frame  stability. 

c.  Mu  may  be  determined  using  linear-elastic 
analysis  (first-order)  combined  with  approximate 
amplification  factors . 
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The  LRFD  Specification  therefore  requires  explicit 
consideration  of  second-order  effects  and  permits  either 
exact  or  approximate  methods  of  analysis . 

Typical  current  design  practice  uses  elastic  first- 
order  analysis  with  approximate  amplification  factors  to 
determine  Mu.  This  procedure  may  be  tedious  and  time 
consuming,  even  if  computers  are  being  used.  In 
addition,  the  approximate  amplification  factors  provided 
in  the  LRFD  Specification,  are  generally  applicable  only 
to  simple,  rectangular  framing.  The  approximate  methods 
of  determining  the  amplification  factors  for  complex 
irregular  framing  are  not  yet  well  developed. 

It  is  the  specific  purpose  of  this  report  to  develop 
an  automated  elastic-second  order  analysis  computer 
program  that  will  meet  the  requirements  of  Section  2a, 
Chapter  H  and  Chapter  C  of  the  LRFD  Manual.  The  use  of  an 
automated  elastic  second-order  analysis  program  will 
remove  the  burden  of  hand  calculations  associated  with 
the  approximate  methods.  It  will  also  permit  accurate 
elastic  second-order  analysis  of  planar  frames  of 
arbitrary  geometry. 


1.3  SCOPE 

1.3-1  Methods  of  Analysis 
The  elastic  second-order  analysis  of  a  beam-column 
originated  from  the  "Buckling  of  Columns"  which  was 
first  introduced  in  1729  by  Van  Musschenbroek  [  2  ]  .  In 
1744,  Euler  introduced  the  first  accurate  stability 
analysis  of  a  centrally  loaded  elastic  column  (  3  ] .  Other 
members  that  were  the  subject  of  subsequent  study 
included  eccentrically  loaded  columns,  columns  with  end 
moments,  and  laterally  loaded  columns.  These  types  of 
columns  were  treated  as  "beam-columns". 

There  are  two  basic  methods  commonly  used  in 
performing  computerized  elastic  second-order  analysis  of 
beam-columns  and  frames .  These  are  called  the  stability 
function  approach  and  the  geometric  stiffness  approach. 
The  latter  is  an  approximate  approach,  and  will  provide 
accurate  solutions  for  a  limited  class  of  problems .  The 
stability  function  approach,  on  the  other  hand,  is 
directly  based  on  the  governing  differential  equations  of 
the  elastic  beam-column,  and  provides  exact  solutions. 
Further,  in  order  to  perform  accurate  elastic  stability 
analysis  requires  the  use  of  stability  functions .  This 
approach  is  used  in  the  development  of  the  Elastic 
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Second-order  Computer  Analysis  Program  [  E  S  C  A  P  ] 
described  in  this  report. 

Timoshenko  and  Gere1  4  ],  Chen  and  Atsuta  [  2  ],  Chen 
and  Lui  [  5  ],  provide  solutions  to  various  cases  of 
elastic  beam-columns .  Salmon  and  Johnson  (  6  ]  give  an 
extensive  list  of  references  on  the  subject  of  beam- 
columns  . 

1.3-2  Assumptions 

The  development  of  -elastic  second-order  analysis  is  based 

on  the  following  assumptions  l  2'  3'  7  3  : 

a.  The  material  is  linearly  elastic. 

b.  Rotations  are  small,  so  that  the  following 
approximations  can  be  made  : 

i.  sin  G  =  6 
ii .     cos  6  =  1 
iii .    62  =  0 
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c.  The  exact  expression  for  the  curvature  of  a  line 
is  : 


1            v" 
4>  =  -  =  *—, T  EQ.      (1) 

P    [l  +  (y')2]3/2 
where  p  =  radius  of    curvature 
<)>  =  curvature 
y  -  deflection 


As  a  consequence  of  the  small  rotation  assumption  (  y  ' 
-  rotation  ) ,  the  usual  relationship  between  curvature 
and  deflection  results  : 

<t>  =  y"  Eq.  (la) 


d.  The  beam-column  is  prismatic. 

e.  Shear  deformation  is  neglected. 

f .  The  frame  members  are  two-dimensional  and 
subjected  to  in-plane  forces  and  deformations 
only. 
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1.3-3   Organization  of  the  Report 

The  report  consists  of  six  chapters.  The  first  chapter 
is  the  Introduction.  Chapter  2  deals  with  the  basic 
theory  of  beam-columns.  The  concept  of  elastic  second- 
order  is  explored  at  some  depth  to  give  an  appreciation 
of  the  subject.  Its  relative  importance  to  the  elastic 
stability  problem  is  briefly  discussed  and  later  covered 
in  Chapter  4. 

Chapter  3  presents  the  methodology  that  is 
implemented  in  the  development  of  the  computer  program  [ 
E  S  C  A  P  ]  .  It  explains  how  an  ordinary  first-order 
linear-elastic  analysis  program  using  the  direct- 
stiffness  method  [  8  ]  is  modified  to  include  second-order 
effects  through  the  implementation  of  stability 
functions.  The  derivation  of  the  beam-column  stiffness 
matrix  is  presented  using  second-order  differential 
equations.  In  Section  3.4-2,  a  method  for  determining 
the  maximum  elastic  second-order  moment  within  a  member 
is  presented. 

In  Chapter  4,  the  reader  is  introduced  to  the 
concept  of  elastic  stability  analysis  of  beam-columns  and 
frames .  The  use  of  Cholesky '  s  decomposition  of  the 
stiffness  matrix  in  determining  the  elastic  buckling  load 
of  frames  and  beam-columns  is  introduced  and  discussed  in 


detail . 

Chapter  5  presents  several  examples  of  beam-columns 
and  frames .  The  computer  solutions  are  compared  to 
closed- form  or  published  solutions .  Some  of  the  computer 
solutions  to  the  sample  problems  are  listed  in  Appendix 
B.  Also,  a  Users  Manual  for  the  computer  program  is 
listed  in  Appendix  A. 

And  finally,  Chapter  6  provides  a  conclusion  on  the 
results  of  the  analysis  of  the  computer  program 

[  E  S  C  A  P  ] .  The  importance  of  elastic  second-order 
analysis  and  elastic  stability  analysis  of  frames  are 
also  summarized.  A  list  of  selected  references  is  also 
provided . 


BACKGROUND  ON  BEAM- 
COLUMNS 


The  term  beam-column  refers  to  a  structural  member 
that  is  subjected  to  both  axial  compression  (or 
tension)  force  and  bending  moment.  It  is  a  structural 
element  simplifying  to  a  beam  or  a  column  as  special 
cases .  A  column  is  a  structural  element  where  axial 
compression  predominates,  while  a  beam  is  a  structural 
element  where  bending  predominates .  The  interaction  of 
the  axial  force  and  bending  moment  may  result  to  a 
nonlinear  behavior  of  a  member  in  terms  of  : 

a.  the  axial  force  and  the  resulting  second-order 
deflections,  and  ; 

b.  the  axial  force  and  the  resulting  second-order 
moments . 

The  elastic  analysis  of  beam-columns  must  consider 
the  effect  of  change  in  geometry  of  the  structural 
element,  i.e,  second-order  effects.  As  shown  in  Fig.  1, 
the  axial  force  will  magnify  the  deflection  and  moment  in 
the  beam.  When  the  axial  force  approaches  a  critical 
value,  the  deflection  and  moment  of  the  beam- 
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(a)    Deflection  due  to  Q 


2nd-order  deflection 

P 


first-order  deflection 


(b)    Deflection  amplified 
by  axial  force  P 


2nd-order  moment 


f  Irst-order  moment 
Mo 


(c)  Moment  amplified 
by  axial  force  P 


Fig.  1     Effects  of  Axial  Force 
in  a  Beam-Column 
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column  increase  without  bound.  This  critical  value  is 
the  elastic  buckling  load  and  can  be  determined  by 
performing  elastic  stability  analysis  of  the  member  or 
structure . 

Most  structures  are  composed  of  many  structural 
elements,  where  each  element  serves  a  unique  purpose,  and 
the  combination  meets  the  functional  objectives  of  the 
structure.  Because  these  elements  are  interconnected, 
they  interact  under  a  given  loading  pattern.  When  both 
axial  forces  and  moments  are  present  in  the  members,  they 
respond  as  beam-columns  and  the  deformations  resulting 
from  the  bending  moments  induced  by  the  externally 
applied  loads  may  be  amplified  by  the  presence  of 
compressive  axial  loads .  However,  when  the  axial  force 
is  tension,  the  effect  on  the  resulting  deformations  is 
opposite  of  the  compressive  axial  force. 

Since  nearly  all  the  members  of  a  structure  are 
subjected  to  both  bending  moments  and  axial  forces,  the 
interactive  effect  of  these  forces  may  be  significant  in 
the  analysis  of  beam- columns .  When  the  magnitude  of  the 
internal  stresses  and  deflections  produced  by  the  axial 
force  in  a  member  is  relatively  small  compared  to  the 
bending  effect,  it  may  be  practical  to  analyze  and  design 
the  member  as  a  beam  element.   On  the  other  hand,  if  the 
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bending  effect  is  relatively  small,  analysis  and  design 
of  the  member  as  a  column  is  appropriate. 

The  presence  of  a  compressive  axial  force  reduces 
the  force  needed  to  cause  a  unit  rotation  or  a  unit 
translation  in  the  transverse  direction  at  one  end  of  a 
member.  Conversely,  the  required  forces  increase  when 
the  axial  force  is  tensile.  In  other  words,  the  presence 
of  axial  force  affects  the  flexural  stiffness  of  a  member 
or  frame.  The  case  in  which  the  axial  force  is 
compressive  is  usually  of  greater  interest,  since  it 
leads  to  the  possibility  of  member  or  frame  instability. 
Using  second-order  analysis  techniques,  the  instability 
can  be  predicted  as  the  axial  compressive  force 
approaches  a  critical  value,  which  reduces  the  flexural 
stiffness  to  zero.  Second-order  analysis  is  therefore 
central  to  the  solution  of  stability  problems. 


3   ELASTIC    SECOND-ORDER 
ANALYSIS 


3.1    GENERAL 

Elastic  second-order  analysis  of  members  and  frames  can 
be  accomplished  by  developing  and  solving  the  second- 
order  differential  equations  of  equilibrium.  Classical 
closed  form  solutions  for  single  members  and  simple 
frames  are  available  l2,4,5,Sl.  However,  for  analysis  of 
large  frames  of  arbitrary  geometry,  closed-form  solutions 
are  difficult  and  cumbersome  to  obtain.  However,  with  the 
advent  of  computers,  automated  second-order  analysis  for 
such  cases  is  possible. 

The  direct  stiffness  method  is  commonly  used  for 
computerized  first-order  structural  analysis .  In  this 
Chapter,  this  method  will  be  used  in  extending  the  first- 
order  analysis  to  a  second-order  analysis  described  in 
the  computer  program  [  E  S  C  A  P  ]  using  stability 
functions . 
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3.2  METHODOLOGY 

As  previously  noted,  the  presence  of  large  axial  force  ( 
i.e,  P  /  P  „  is  large,  where  P  .  is  the  elastic  buckling 
load  )  /  affects  the  geometry  and  stiffness  of  a  beam- 
column  element.  If  the  axial  force  P  is  known,  the 
maximum  deflection  and  bending  moment  can  be  calculated 
using  second-order  formulations .  However,  the  member  or 
frame  is  initially  indeterminate  and  the  axial  forces  are 
unknown.  The  axial  forces  are  related  to  the  joint 
displacements  and  must  be  first  calculated  because  the 
axial  forces  influence  the  members'  fixed-end  actions, 
which  are  the  elements  of  the  load  vector  used  in 
determining  the  joint  displacements.  Additionally,  the 
axial  forces  also  influence  the  stiffness  matrix  of  the 
member  or  frame.  One  way  of  determining  the  initial 
axial  forces  is  to  perform  a  first-order  linear-elastic 
analysis  (  first  cycle  )  using  the  direct-stiffness 
method.  This  will  give  the  first-order  shears,  moments, 
deflections  and  axial  forces . 

In  the  second  cycle,  the  previously  computed  axial 
forces  are  used  to  modify  member  stiffness  and  fixed-end 
actions .   Then  new  values  of  joint  deflections  and 
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Fig.  2  -  Flowchart 
Elastic  Second-Order  Analysis 
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member-end  actions  are  calculated.  This  process  is 
repeated  until  the  differences  between  the  deflections  in 
successive  cycles  are  within  the  acceptable  tolerance. 
This  is  the  basic  algorithm  used  in  implementing  the 
computer  program  [  E  S  C  A  P  ] .  The  flowchart  shown  in 
Fig.  2  outlines  the  basic  algorithm  of  the  computer 
program . 


3.3   DERIVATION  OF  THE  SECOND-ORDER  STIFFNESS  MATRIX 

3.3-1  Translation  Stability  Functions 
To  begin  our  derivation  of  the  stiffness  matrix  of  a 
beam-column,  let  us  consider  Fig.  3  .  The  beam-column 
element  is  subjected  to  a  unit  translation  (  Fig.  3a  ) 
and  a  unit  rotation  (  Fig.  3b  )  under  the  presence  of  a 
compressive  axial  load  P  The  elastic  stiffness 
coefficients  for  these  two  beam-column  elements  will  be 
derived  from  second-order  differential  equations  of 
equilibrium. 

The  degrees  of  freedom  of  a  beam-column  element  are 
shown  in  Fig.  4a.  A  free-body  diagram  of  a  section  of 
the  beam-column  element  cut  at  distance  x  from  the  left 
support  is  shown  in  Fig.  4b.   The  sum  of  external 
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(a)  unit  translation 


(b)  unit  rotation 


Fig.  3       Stiffness  Coefficients 
of  a  Beam-Column 
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(a)    member  degrees  of  freedom 


(b)    free-body  diagram 
unit  translation 


(c)  free-body  diagram 
unit  rotation 


Fig.  4       Unit  Translation  -  Unit  Rotation 
of  a  Beam-Column 
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moments  at  point  0  (Fig.  4b),  assuming  counter-clockwise 
moment  is  positive  is  given  by 

Moxt  -  -P(A  -  y)    +  SM22  -  SM22x  ,  Eg.    (2) 

where   SM  i3   =  is   the  member  stiffness   coefficient 
corresponding  to  the  degree  of  freedom  i   induced  by  a 
unit  displacement  corresponding  to  the  degree  of  freedom 
j    while  the  other  displacements  are  kept  to  zero. 
Applying  the  differential  equation  of  an  elastic  curve 

M  =  -Ely"  ,  Eg.  (3) 

then  the  interior  moment  must  balance  the  exterior 
moments,  resulting  in  : 


M  =  -EI^-l.   = 
dxz 

P(A  -  Y)    +  SM32  -  SM22X  Eg.    (4) 

M  =  EI-^  =  P  (A  -  Y)    +  SM22x  -  SM32 
dx2 


Eg.    (5) 
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Rearranging  the  terms  of  Eq.  (  5  )  yields 


°-JL   +  k2y  =  k2A   +  ^ ^-    Eg.    (6 

dx2  ^J" 


p 
where  k2   =  —  Fg.  (7) 

EI 


The  complementary  solution  of  Eq.  (  6  )  is  given  by 


yc  =  A  sinkx  +  B  coskx  Eg.    (8) 


and  the  particular  solution  of  Eq.  (  6  )  has  the  form 

yp  =  Cx  +  D  .  Eg.    (9) 

Differentiating  Eq.  (  9  )  twice  with  respect  to  x  and 
substituting  the  first  and  second  derivatives  into  Eq.  ( 
6  )  yields  Eq.  (  10  )  .  The  constants  C  and  D  are 
determined  by  comparing  coefficients  of  the  terms  on  the 
left  hand  side  of  Eq.  (  10  )  to  the  terms  on  the  right 
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hand  side  as  shown  below 


0  +  k2{Cx  +  D)    = 


.2  A  SM22X     .        SM22 


k2A     +     J^ZZJ_      _     J^3Z  £g  (10) 


where  C 


=      *fc    ;    D  =   A    - 


£2£J  ic2EJ 


The  general  solution  of  Eq.  (  6  )  is  the  sum  of  its 

particular   solution    and   complementary  solution. 

Inserting  the  values  of  C  and  D  into  Eq.  (  9  )  ,  the 
general  solution  has  the  form 


y  =  A  sinicx  +  B  coskx  +  A  + 

Eq.    (11) 


{SM22x  -  SM32) 


k2EI 
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To  solve  for  the  constants  A  and  B,  apply  the  boundary 
conditions  at  joint  k  of  Fig.  4b  .  At  x  =  0,  y  =  A; 
then 


B  =  —3!  .  Eg.    (12) 


k2EI 


At  x  =  0,  dy/dx  =  0;  then 


A  =  -—22-  .  Eg.  (13) 

Je3EJ 


The  stiffness  coefficients  SM  22  and  SM  32  are  determined 
by  applying  the  boundary  conditions  at  joint  m  of  Fig. 
3a.   When  x  =  L  ,  y  =  0  ;  then 


A  sinkL   +  B  coskL   +  A  +  ^— 21_  =  o 

k2El 


Eg.    (14) 


When  x  =  L  ,  dy/dx  =  0  ;  then 


Ak  coskL  -  Bk  slnkL   +  21  =  o      Eg.    (15) 

k2EI 
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Inserting  the  values  of  A  and  B  into  Eqs .  (14)  &  (15) 
yields  Eqs.  (16a)  &  (16b)  : 


SM22 


sin  kL  +  —  cos  kL  +  A  + 


k3EI  k2EI 


{SM22L  -  SM„) 

32-  =  0         Eg.    (16  a) 


k2EI 


SM22  ,     SM,2       .   , 

cos  kL  -  -r^   sin  kL 


k2EI  kEI 


SM*2    =  0  .  EQ.    (16b) 


k2EI 


Simplifying  Eqs.  (16a)  &  (16b)  yields 


SM22  (kL  -   sin  kL)    + 
SMsz  k(cos  kL  -   1)  =  -k2ElA  Eg.  (17) 

SM22   (1  -  cos  kL)   -  SM32   sin  kL  =  0     Eg.    (18) 


Solving  Eqs.  (17)  &  (18)  for  SM  22  yields 


SM22    [kL  sin  kL   -  sin2  kL   + 

(cos  kL  -  1  -  cos2  kL  +  cos  kL)  ] 

=  -k3 EI  sin  kL    A  Eg.  (19) 


Simplifying  Eq.  (19)  yields 


Si*^    [2-2  cos  kL  -  kL  sin  kL]    = 
ir3£J  sin  JcL    A  Eg.    (20) 


24 


Multiply  Eq.  (20)  by  12L  3  /  12L 3  and  introduce  the  term 


<t>  =  (2-2  cos  kL   -  kL  sin  kL)  Eg.    (21) 
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the   stiffness   coefficient   SM 22      is 


SMi2  =   i*S)  {6kL^sinkL)        EQ     (22) 
L3  12  <pc 


SM27.  =  -^4^(6  Si)   A  Eq.    (23) 


where  S  x   is  referred  to  as  a  stability  function  and  is 
equal  to 


=  (JcLpsinJeL     Eq^    (24) 
12  $c 


Substitute  Eq.  (22)  into  Eq.  (17)  and  solve  for  SM  32 


•^32  =  ^tt^  (3L  s2)  A    Eq.    (25) 


The  other  translational  stability  function  S  2  is 


,  (*L)»  (l-cos*L)    £   (26) 
6  ♦<, 
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3.3-2   Rotation  Stability  Functions 

Consider  now  the  beam-column  element  in  Fig.  3b  which  is 
subjected  to  a  compressive  axial  force  P  and  a  unit 
rotation  at  joint  n  .  A  free-body  diagram  is  taken  at 
a  distance  x  from  the  left  support  as  shown  in  Fig.  4c 
Summing  the  external  moments  about  point  0,  gives 


Mm*   =  Py  -  SM23*   +  SAf33    Eg.   (27) 


Applying  the  differential  equation  of  an  elastic  curve 
{Eq.  (2)  }  and  using  equilibrium,  the  external  moments 
{Eq.  (27)},  must  be  equal  to  the  internal  moment. 
Rearranging  the  terms  yields 


EI^-^-   +  Py  =  SM22x  -  SM22  Eq.   (28) 

dx2 
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The  complementary  solution  of  Eq.  (28)  is  the  same  as  Eq. 
(8)  .  Using  similar  procedures  described  in  Section  3.2- 
1,  the  particular  solution  of  Eq.  (28)  has  the  form 


yB  =  -r-^-x   -  21     Eq.    29 

y*        k2EI  k2EI 


The  general  solution  of  Eq.  (28)  is  the  sum  of  the 
complementary  solution  {Eq.  (8)  }  and  the  particular 
solution  {Eq.  (29) )  ,  and  has  the  form 


y  =  A  sin  kx  +  B  cos  kx  +  —  x 

kzEI 


5^3  Eq.  (30) 


k2EI 


„/_*]_    —  ^0      T,v    _     DL.    »in      Irv    J.  "23 


y'  =  Ak  cos  kx  -  Bk  sin  ice  +  —        £g.    (31) 

k2EI 
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The  constants  A  and  B  of  Eq.  (30)  are  determined  by 
applying  the  boundary  conditions  at  joint  n  of  Fig.  3b: 


OJLf 

at  x  =  0  ;  y  =  0  then    B   =  —  Eg.    (32) 

k2EI 


fl    SM 

at  x  0  ;  y'  =  6  then    A   =  -f  -  ^-         Eq.    (33) 

k        k2EI 


Inserting  the  values  of  A  and  B  to  Eq.  (31)  yields 


y  -  (-r  _  — — )  £  cos  iex  -  —  k  sin  kx  + 


-^       Eg.  (34) 
£2EJ 
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Applying  the  boundary  conditions  (  joint  r  of  Fig.  3b 
)  to  Eq.  (  34  )  at  x  -  L  ,  y'  =  0  ;  then 


(  6  _  _S^L)   k  CQS  kL  _    SM^L  k  sin  kL   + 
k       k2EI  k2EI 


— -22.)  =  0      Eq.    35 
k2EI 


At  x  =  L  ,  y  —   0  ;  then 


(  D   _  ^3_j   s±n  kL    +  33_  CQS  ^  + 

ic   Jc3£-J  k2EI 


<J"k)  L-  J^l=0       *g.(36) 
fc2FJ      Jc2£J 
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Solving  Eqs.  (35)  and  (36)  simultaneously  yields 


Elk 

4>c 


5Af33  =  -==£  [sin  kL  -  kL  cos  kL]    *  6 


Eg.    (37) 
where    <J>_  =   [2-2  cos  icL  -  kL  ]     Eg.    (21) 


Multiplying    Eq.     (37)    by    4L3    /    4L3    and    simplifying    the 
terms    : 


SAT      =  Ml   *    (2L2)    *    (s3)   0  Eg.  (38) 

L3 

wAere 

icL   (sin  kL  -  kL  cos  icL)  „„     ,-,QX 

s3  -  — —t r,g.    ijy; 

4  <j)c 


Substitute  Eq.  (37)  into  Eq.  (36)  yields 


SM23   =  -^-^-  (  1  -  cos  kL   )  0     Eg.  (40) 
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Equations  (25)  and  (40)  are  equivalent  and  satisfy 
Betti's  Law  and  Maxwell's  Law  (  X1  ] .  This  indicates  that 
the  second-order  stiffness  matrix  is  symmetric.  The 
remaining  rotation  stability  function  is  calculated  by 
summing  moments  at  joint  r  of  Fig.  3b  assuming  counter- 
clockwise moment  is  positive.   By  equilibrium  we  have 


SM23  L  -   P(0)  +  SM23   +   5A^3  =  0      Eg.   (41) 


Substituting  Eqs.  (37)  &  (40)  into  Eq.  (41)  yields 


SM63  =    EI£2Q  [L  -  L  coskL]    - 


^P^  [sin  kL  -  kL  cos  kL] 
-^^    [kL  -  sin  kL  ]  Eg.  (42) 
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Introducing  the  term  2L  3/2L3  in  Eq.  (42)  and  simplify 
to  Eq.  (43) ,  which  is  the  coefficient  of  stiffness  at 
joint  r  caused  by  a  unit  rotation  at  joint  n  ,  Fig.  3b  : 


5^3  =  ^r   *   2I.2  (s4)  *  6     Eq.  (43) 


where 


=  kL   jkL   -  sin  kL)       EQm    (44) 
2  $c 


For  negative  values  of  axial  force  P  (tension)  , 
the  stability  functions  are  derived  in  similar  fashion 
applying  the  same  concept  of  elementary  beam  analysis .  To 
illustrate,  let  us  investigate  Eq.  (4)  .  Assuming  that 
axial  tensile  force  is  negative,  then  the  homogeneous 
second-order  differential  equation  has  the  form  of  Eqs . 
(45)  &  (46) 
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Ely"  -  Py  =  0  Eg.    (45) 

y"  -  k2y  =  0  Eg.    (46) 

where  k2   =  —     Eg.    (7 ) 
EI 


The  complimentary  solution  of  Eq.  (46)  is 

yc  =  A  sinh  lex  +  B  cosh  iex     E<y.  (47) 
yc  -  A  sin  lex  +  B  cos  lex     .Eg.  (8) 


Equations  (8)  and  (47)  are  similar  in  form  except  the 
former  contains  the  trigonometric  functions  while  the 
latter  has  the  hyperbolic  functions.  Using  the 
complimentary  solution  of  Eq.  (47)  ,  the  stability 
functions  for  an  axial  force  in  tension  are  calculated 
following  the  same  procedures  as  when  the  axial  load  is 
in  compression.  The  four  stability  functions  are  listed 
in  TABLE  I  for  axial  forces  in  compression,  tension,  and 
zero  I12] . 


TABLE  -  I   STABILITY  FUNCTIONS  34 

OF  A  BEAM-COLUMN 


Stability   Functions   for  a  Beam  Subjected  to 
Axial  Forces 

Direction  of  Axial  Force 

Function   Compression  Tension 

s  i        (kL) 3  sin  kL  (kL) 3  sinh  kL 


12  * c  12  ^  t 

(kL)2  (1  -  cos  kL)      (kL)2  (cosh  kl  -  1) 


6  <J>  c  6  4  t 


s  3    kL  (sin  kL  -  kL  (sinh  kL  - 

kl  cos  kL  kl  cosh  kL) 


4  <>e  4  +t 

s  4    kL  (kL  -  sin  kL)         kL  (sinh  kL  -  kL) 


2  4>c  2  <K 

<}>  e  =2-2  cos  kL  -  kL  sin  kL 
<t>  t   =2-2  cosh  kL  +  kL  sinh  kL 
Note:  When  axial  force  is  zero  : 
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The  member  stiffness  matrix  with  the  stability 

functions  is  represented  by  Eq.  (48),  where  [  KB  ]  is  the 

modified  stiffness  matrix  and  includes  the  effects  of 

axial  forces  in  tension  or  compression. 


uy  -  4 


A 

0 

0 

-A 

0 

0 

0 

12X 
L2     1 

^s2 

L     2 

0 

"  12  J*x 

L2     l 

^s2 

L     2 

E 

0 

L*2 

4Js3 

0 

~^s2 
L     2 

2 1  S< 

L 

-A 

0 

0 

A 

0 

0 

0 

12T 

'  l2^ 

~^s2 

L     2 

0 

12Is 

~^s2 

L     2 

0 

^s2 

L     2 

21s4 
Eq. 

0 
(48) 

-^s2 

L     2 

4Js3 
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When  the  axial  force  is  zero,  the  stability  functions  are 
reduced  to  unity  and  Eq.  (48)  is  reduced  to  Eq.  (49) 


[*.]  =  T 


A 

0 

0 

-A 

0 

0 

12J 
L2 

6J 
L 

0 

12J 
L2 

0 

6J 
L 

4J 

0 

_6T 
L 

-A 

0 

0 

A 

0 

0 

12J 
L2 

_6J 
L 

0 

12J 
L2 

0 

6J 
L 

21 
EQ. 

0 
(49) 

_6J 
L 

0 

6J 
L 

21 

0 

_6J 
L 

4J 


Equation   (49)   is  the  first-order  elastic  stiffness 
matrix . 


3.4   LOAD  FUNCTIONS  AND  MAXIMUM  MOMENTS 


3.4  -  1 


Load  Functions 


The  fixed-end  actions  for  a  beam-column  produced  by 
applied  lateral  loads  depends  upon  the  magnitude  of  the 
axial  force  acting  on  the  member  as  well  as  upon  the 
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distribution  and  intensity  of  the  lateral  loads .  These 
fixed-end  actions  or  joint  loads  are  used  in  determining 
the  joint  deflections  and  member  forces  in  the  structure . 
The  fixed-end  actions  form  the  elements  of  the  load 
vector  in  the  direct  stiffness  method. 

The  two  types  of  member  loading  (  Fig.  5  ) 
implemented  in  the  computer  program  are:  (1)  uniform  load 
distributed  over  the  full  length  of  the  member,  and  (2) 
concentrated  loads  applied  anywhere  along  the  span.  The 
procedures  used  in  deriving  these  fixed-end  actions  are 
similar  to  those  described  in  Section  3 . 3  and  will  not  be 
repeated  here.  A  complete  derivation  of  the  fixed-end 
actions  for  the  two  types  of  loading  is  presented  in 
detail  in  reference  [  14  ]  .  Timoshenko  and  Gere  [  4  ] 
present  solutions  to  specific  cases  of  beam-columns  with 
various  constraints  at  the  ends .  This  includes  a  member 
simply  supported  on  one  end  and  built  in  at  the  other 
end,  as  shown  in  Example  5.2-3,  Section  5.2  .  Ghali  and 
Neville  [  9  ]  present  a  different  procedure  in  determining 
the  fixed-end  actions  of  a  beam-column  using  matrix 
methods . 

The  fixed-end  actions  of  a  beam-column  for  the  two 
types  of  member  loading  mentioned  above  are  shown  in 
figures  Fig.  6  through  8  [  10'  12'  14  ]. 
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(a)  Concentrated  Load 
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(b)  Uniform  Load 


Fig.  5         Types  of  Loading 
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Axial  Compression 


FEA^  ■  -j-   (cosier,  -  coska  +  cosicb  +  kb  sinJeL  -  1) 


FEA,  -  -r?-  ( sin  JcL  -  sin  ka  -  sin  kb 
-  kb  cos  JcL  +  kL  cos  Jcb  -  ka) 


FEAs  *  -^  (cos  JcL  -  cos  Jcb  +  cos  Jca 
+  Jca  sin  JcL  -  1 ) 


FEA*  ■  — =?-  ( -sin  JcL  ♦  sin  Jcb  +  sin  Jca 

+  Jca  cos  JcL  -  JcL  cos  Jca  +  Jcb) 
$c  s  2  -  2  cos  JcL  -  JcL  sin  JcL 


FEA 


Concentrated  Load 
Fig.    6  -  Fixed-End.  Actions  for  a  Concentrated  Load 
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Axial  Tension 


=  0_ 
4>t 


FEA2   =  -p-  ( cosh  kL  -   cosh  ka  +   cosh  kb 


-  kb  sinh  kL  -  1 ) 


ic4>t 


FEA3  =  -j-^-    (sinh  kL  •-  sinh  ica  ■-  sinh  kb 


-  kb  cosh  kL  +  kL  cosh  kb  -  ka) 
FEAr  =  -r-  (cosh  kL  -  cosh  icb  +  cosh  ka 

5      4>t 

-  ka  sinh  ieL  -  1 ) 

FEA,  =  -ri-  (  -sinh  icL  +  sinh  kb  +  sinh  ica 
^        ic<|)t 

+  ica  cosh  kL  -  kL  cosh  ica  +  icb) 
4>t  =  2  -  2  cosh  icL  +  icL  sinh  icL 


*  =  N  — 
M  £ J 


Fig.   7  -  Fixed-End  Actions  for  a  Concentrated  Load 
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Axial   Compression 


FEAj 


FEflj  =  FEAg 


wL2    r      12       + 


2 

6  sinJcL 


12         (JcL)2      JcL  (cos  JcL  -  1) 


] 


FEA^ 


Axial  Tension 


FEflj  -  Ffiflg  «  - 


wL2  r      12 


wL 


6  sinh  JcL 


12        (JcL)2        *L  (cosh  JcL  -  1) 


FEflg  ■  -J!Efl3 


FEA 


2      l    ^r 


Uniform  Load 


Fig.  8  -  Fixed-End  Actions  for  a  Uniform  Load 
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3.4-2  Maximum  Elastic  Second-Order  Moments 

It  is  important  for  a  designer  to  know  the  location  and 
magnitude  of  the  maximum  moment  within  a  member.  Some 
structural  analysis  programs  only  calculate  the  member 
end  actions  (  i.e.  shears, moments  and  axial  forces  )  and 
the  joint  deflections.  In  some  cases,  the  maximum 
moments  occur  at  the  joints,  when  the  axial  forces  are 
small  and  no  member  lateral  loads  are  present.  However, 
when  the  axial  forces  acting  on  the  members  become  large, 
the  maximum  moment  may  occur  within  the  span  as  shown  in 
Fig.  9  ,  even  if  no  lateral  loads  are  present. 

In  this  section,  the  derivation  of  the  maximum 
moment  is  limited  to  the  case  of  a  uniform  distributed 
load.  The  formulation  for  this  case  is  straight  forward, 
since  an  equation  for  bending  moment  can  be  expressed  as 
a  continuous  function  along  the  length  of  the  member. 
The  bending  moment  at  a  concentrated  load  point  can  be 
easily  obtained  by  introducing  a  node  at  the  location  of 
the  load. 

The  maximum  moment  in  a  member  with  a  uniform 
distributed  load  is  calculated  based  on  the  procedure 
suggested  by  White  and  Ha j jar  [  13  ]  .   An  isolated 
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Fig.  9       Maximum  Elastic 
Second-Order  Moment 
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member  of  a  frame  structure  is  treated  as  a  simply 
supported  beam-column  that  has  joint  rotation  and  no 
translation . 

To  illustrate  this  procedure,  Fig.  10  shows  a  rigid 
frame  with  its  deformed  shape  when  subjected  to  the  two 
types  of  loading  described  in  Section  3.4-1.  Member  1 
is  isolated  and  the  free-body  diagram  is  shown  in  Fig. 
10b  . 

Let  us  assume  the  member  end  actions  are  known  from 
the  results  of  the  elastic  second-order  analysis .  When 
we  connect  joints  1  and  2  <  Fig.  11a  ),  the 
member ' s  chord  location  is  established  and  calculated 
using  the  new  joint  coordinates  If  the  member  end 
forces  and  any  other  loads  applied  directly  to  the  member 
are  transformed  to  local  axes  oriented  along  the  member '  s 
chord,  White  and  Ha j jar  suggest  that  there  is  no 
difference  between  the  behavior  of  this  isolated  member 
(  Fig.  10  and  11  )  and  the  behavior  of  a  simply  supported 
beam-column  with  the  same  applied  loads  (  Fig.  15, 
Example  5.2-1  ). 

To  derive  the  location  of  the  maximum  elastic 
second- order  moment,   we  shall  use  the  differential 
equation  of  an  elastic  curve  {  Eq.   (3)   }  about  the 
transformed  local  axes  of  Fig.  lib  . 
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Loading  and 
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Fig.  10         Derivation  of  Maximum 
Elastic  Second-Order  Moment 
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V    =  V  cos  +  P'  sin 

a 

(a)  Element  end  forces 

oriented  at  the  chord's  axis 


Ma 


t— — I 


a 


(b)  free-body  diagram 


Fig.  1 1        Transformed  Axis 
of  a  Beam-Column 
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The  second-order  differential  equation  of  a  beam-column 
subjected  to  a  uniformly  distributed  load  with  member  end 
forces  known  is  given  by 


Bigz  -  -  ut  -  P(-y)  -  vmx  *  Jc*L*l: 

dx2  2 


Eg.    (50) 


The  general  solution  of  Eq.  (50)  is  given  by 


y  -  A  siniex  +  B  coskx  +  — = —  [ 

k2EI        2 


*  Vax  -  JL  -  Ma]  Eg.    (51) 

k* 
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Differentiating  Eq.  (51)  twice  and  applying  the 
differential  equation  of  an  elastic  curve  {  Eq.  (3)  } 
yields 


y"  -  -A  k2   sin  kx  -  B  k2   cos  kx 


+ 


w 


k2EI 
Eq.    (52) 

Mx   =  -Ely" 

-  Elk2  [A  sin  kx  +  B  cos  kx]    -   — 

k2 

EQ.   (53) 


When  x  =  0 ;  Mx  =  Ma 

then  B  =  \[Ma  +  —}  .  Eq.    (54) 

P  k2 

When  x  =  L  /  Mx  =  Mb 

Mh  -  M.  cos kL  w  jcL    „   ,„i 

then  A  =  -^ ?__; +  — —  tan4r   Eq.    (55) 

PsmkL  pk2  2 


where    P  =  k2EI 
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Differentiating  Eq.  (53)  once  with  respect  to  x  gives 
the  expression  of  shear  force  anywhere  along  the  span  and 
has  the  form 


VL  =  m'x  =  [A  coskx  -  B  sinicx]  Elk2 


Eg.    (56) 


When  the  shear  force  is  equal  to  zero,  the  moment  is  at 
local  maximum.  Inserting  the  values  of  A  and  B  into  Eq. 
(56) ,  the  location  of  the  maximum  moment  anywhere  along 
the  span  of  a  beam-column  with  a  uniform  distributed  is 
determined  using  Eq.  (57) : 


tan  kx  = 

w            kL      Mh  -  AT.  cosieL 
k2            2                sinieL 

MM  +  — 
a         k2 

Eg.    (57) 


where  x  is  measured  from  the  left  support.  All  the 
applied  loads  must  be  transformed  in  the  local  axes  when 
using  Eq.  (57) . 

When  the  end  moments  in  Fig.  11,  are  equal  to  zero, 
the  constants  A  and  B  are  given  by  Eq.  (58)  .  The  moment 
anywhere  along  the  member  is  then  given  by  Eq.  (59) 
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B  =  -^-     ;  A  =  — ^tan  -^      Eg.    (58) 
Pk2  Pk2  2 


Mx  =  -^  [tan  -^  sin  icx  +  cos  Joe  -  1] 
k2  2 


£<J.  (59) 


The  maximum  moment  for  a  simply-supported  beam-column 
with  uniform  distributed  load  is  located  at  the  midspan. 
Thus  , 


kL  L 

tan  kx  -   tan —  ;  where    x  -  —   .  Eg.    (60) 

2  2 


Inserting  the  value  of  x  into  Eq.  (59)  yields 


Mx  =  —    [tan  i^sin  ^  +  cos  ^  -1] 
x        k2  2  2  2 


Eg.    (61) 
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Introducing  the  term  u  ,  Eq.  (61)  becomes 


Mx  =  —  [tanu  sinu  +  cosu  -  1]     Eg.    (62) 
k2 


where    u  =  — ^= 


Applying  trigonometric  identities  and  simplifying  the 
terms  in  Eq.  (62)  yields 


Mx=      »£!  [  2  (sec  u  -   X)  ]  .     Eg.    (63) 
8        u2 


Equation  (63)  is  the  general  expression  of  the  maximum 
moment  at  midspan  of  a  simply-supported  beam-column  with 
uniform  distributed  load.  The  first-order  moment  (  wL 2 
/  8  )  is  amplified  by  the  terms  inside  the  brackets . 
This  term  is  called  the  moment  amplification  factor 
(MAF)  . 

Let  us  now  investigate  what  happens  when  the  lateral 
load  w  =  0  and  M  ,  ,  M  b  are  not  equal  to  zero . 
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Applying  the  boundary  conditions  at  x  =  0  and  x  =  L,  from 
Eq.  (51)  ,  the  constants  A  and  B  are  determined.  When  x 
=  0  ,  y  =  0  ,  and  w  =  0  ;  then 


B  =  -^  .  Eg.    (64) 


When  x  =  L  ,  y  =  0  ,  and  w  =  0  ;  then 


A  *  fh  tan  ^  -    „».    Z\  Eg.    (65) 

P  2         PsinicL 


where  VBL  =        a       b  Eg.    (66) 


Inserting  the  value  of  the  trigonometric  identity  {  Eq. 
(67)  }  into  Eq.  (65), 


tan  *£  =  (1  -.cosJcL)  (67) 

2       sin£L 
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Then  the  constant  A  is  given  by  Eq.  (68) . 


{M.COSkL  +  AL) 

A=~ »     •     ,  ,  EQ.    (68) 

PsmkL 


Inserting  the  values   of  A  and  B  into  Eq.     (53)    when     w  is 
equal      0    ,    yields 


A  (Mh  +  Af,cos  kL) 

tan  kx  =  ^  =  -  — b        a        tfg.    (69) 

S  AfasmicL 


«_      =  -    [ — 2 _J_^ L]  sinicx  +  Mmcoskx 

x  smicL  a 

Eg.  (7  0) 


Thus,  Eq.  (70)  is  the  moment  at  any  location  along  a 
beam-column  when  subjected  to  end  moments  M  ,,  M  b  and  an 
axial  force  P.  Inserting  the  value  of  kx  from  Eq.  (69) 
into  Eq.  (70) ,  the  maximum  moment  is  given  by  Eq.  (71) . 
The  minus  sign  that  appears  in  Eq.  (71)  simply  indicates 
that  the  moment  M  max  causes  tension  on  the  top  fiber  of 
the  cross  section  [  5 ] . 
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^nax    =    ~Mt 


(Ma/Mb)2  +  2    {Mm/Mb)  coskL  +  1 
N  sin2  JeL 

£g.    (71) 


When  end  moments  are  present  in  the  member,  there  is 
always  a  possibility  that  the  maximum  moment  will  occur 
at  the  end  of  the  span.  Therefore,  the  value  of  x  in  Eq. 
(69)  must  always  be  evaluated.  When  x  is  greater  than 
the  span  length  or  less  than  zero,  the  maximum  moment 
occurs  at  the  end  of  the  member  and  it  is  equal  to  the 
larger  of  the  two  end  moments . 

Equations  (51),  (53),  (54),  (55)  and  (57)  are  used 
in  the  computer  program  to  calculate  the  location  of  the 
maximum  moment  of  a  member  with  uniform  load  and 
subjected  to  end-moments .  These  equations  have  been 
developed  for  the  case  where  the  axial  force  is  in 
compression.  For  axial  force  in  tension,  similar 
equations  are  used,  except  hyperbolic  functions  replace 
the  corresponding  trigonometric  functions.  The  maximum 
moment  is  then  compared  with  the  end  moments  of  the 
member,  and  the  larger  of  the  three  controls.  The 
program  also  calculates  the  deflection  at  the  point  of 
maximum  moment.   If  the  maximum  moment  occurs  at  one 
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end  of  the  member,  then  the  computer  program  prints  out 
the  deflection  at  that  joint.  For  deflections  within  a 
member,  the  computer  program  only  calculates  the 
displacements  relative  to  the  chord  member.  (See  Example 
5.2-3) 


ELASTIC   BUCKLING   LOAD 


4  . 1       GENERAL 

This  chapter  deals  with  the  elastic  stability  analysis  of 
beam-columns  and  frames .  The  same  stability  functions 
and  stiffness  matrix  developed  in  Chapter  3  are  used  in 
determining  the  elastic  critical  loads.  A  load  factor 
increment  is  applied  successively  to  the  external  loads 
until  the  stiffness  matrix  of  the  structure  has  become 
singular.  The  decomposition  of  the  stiffness  matrix 
using  Cholesky' s  method  will  also  be  presented.  Several 
examples  are  examined  and  presented  in  Chapter  5. 


4.2   COMPUTER  ANALYSIS  OF  THE  ELASTIC  STABILITY  OF 
PLANE  FRAMES 

In  Section  3.3,  the  stability  functions  of  the  beam- 
column  were  derived  and  used  to  modify  the  member 
stiffness  matrix  {Eq.  (48) }  to  include  the  effects  of 
axial  forces.  The  joint  displacements  and  member  end 
forces  are  calculated  iteratively. 
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The  computer  program  was  modified  to  calculate  the 
elastic  critical  load  using  load  increments.  The  elastic 
critical  load  is  determined  by  performing  successive 
computations  until  the  stiffness  matrix  of  the  structure 
becomes  singular. 

In  the  computer  program,  the  structure's  stiffness 
matrix  is  decomposed  using  Cholesky's  method  l  8'  14  ]. 
The  stiffness  matrix  [  KB  ]  {Eq.  (48)}  is  factored  into 
an  upper  and  lower  triangular  matrices . 

In  general,  structural  systems  generate  stiffness 
matrices  that  are  symmetrical,  square,  and  with  positive 
diagonal  elements.  They  display  the  unique  property  of 
positive-definiteness.  If  this  is  the  case,  then 
Cholesky's  method  is  valid  and  the  stiffness  matrix  can 
be  expressed  as  {  K B  }  =  {U  fc}  *  {  U  } .  The  matrix  {  U 
}  represents  the  upper  diagonal  elements  of  the 
structure's  stiffness  matrix  {K  .  }  after  the 
decomposition  process .  The  elements  of  the  lower 
triangular  matrix  are  represented  by  the  matrix  {  D  t  } , 
which  is  the  transpose  of  the  matrix  {  U  }  .  Eq.  (72) 
shows  the  elements  of  the  structure's  stiffness  matrix  { 
K m  }  and  the  elements  of  the  upper  and  lower  triangular 
matrices  after  the  decomposition  process . 
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Eg.    (72) 


Performing  the  indicated  multiplication  of  {  U  }  and  {  U 
t  }  and  equating  the  product  to  the  elements  in  {  K  m  } 
yields  Eqs .  (73)  &  (74)  shown  in  Fig.  12.  If  during 
decomposition  one  of  the  diagonal  elements  u  X1  approaches 
zero  or  is  less  than  zero,  the  structure  has  reached  the 
elastic  buckling  load. 
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where  u„  = 


ii 


N 


««i-l 

*ii  "     £    u^l         e<2-    <73> 


a-l 


are  diagonal  elements  of  [U] 
k±1  =  diagonal  elements  of  [K] 


a  =  i-l 


u 


ij 


a  =  l 


Eg.    (74) 


zii 


i  i  J   (Jo^er  diagonal  elements) 
i  <  j   (upper  diagonal  elements) 


Fig.  -  12   DECOMPOSITION  OF  MATRIX  [  K.  ] 
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To  better  understand  the  concept,  a  flowchart  is 
shown  in  Fig.  13.  In  the  flowchart,  the  process  starts 
with  the  calculation  of  the  first-order  axial  forces. 
Then  these  axial  forces  are  used  to  modify  the 
structure's  stiffness  matrix  [  K B  ] .  The  subroutine 
DECOMP  decomposes  the  structure's  stiffness  matrix  using 
the  procedure  described  in  Fig.  12.  The  program  then 
checks  the  diagonal  elements  u  u  for  value  less  than  or 
equal  to  zero.  If  all  the  diagonal  elements  are  greater 
than  zero,  a  new  load  factor  (  LF  )  is  calculated  and 
the  process  continues . 

If  any  one  of  the  diagonal  elements  is  less  than 
zero,  the  program  stores  the  load  factor  LF,  calculates 
the  joint  displacements,  prints  the  output  and  then 
terminates  the  process .  The  output  contains  the  load 
factor  LF  and  the  joint  displacements .  To  obtain  an 
accurate  buckling  load,  small  load  increments  must  be 
used. 

The  buckling  mode  shape  is  determined  by  converting 
the  calculated  joint  displacements  into  relative  joint 
displacements.  The  joint  displacements  at  buckling  are 
expressed  in  terms  of  relative  displacements . 

This  approach  of  determining  the  elastic  buckling 
load  is  simple  to  implement  in  the  computer  program, 
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Fig.  13  Elastic  Stability 

Analysis  of  a  Frame  and  Beam-Column 
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since  the  decomposition  of  the  matrix  is  routinely 
performed  during  the  first  and  second-order  analyses . 

A  flowchart  showing  the  subroutine  DECOMP  '  8  ]  is 
shown  in  Fig.  14.  The  decomposed  structure  stiffness 
matrix  is  also  used  in  solving  for  the  unknown 
displacements.  In  the  flowchart,  there  are  two  returns. 
One  occurs  when  the  value  of  I  becomes  greater  than  N  and 
the  processing  is  complete.  The  second  is  when  the 
decomposition  fails .  The  latter  serves  as  a  safeguard 
against  an  undefined  arithmetic  operation  and  it  also 
indicates  that  the  critical  load  of  the  structure  has 
been  reached.  The  undefined  arithmetic  operation  will 
result  to  an  error  in  determining  the  critical  load  if  at 
the  beginning  of  the  linear-elastic  (first-order) 
analysis,  the  structure  stiffness  matrix  is  not  positive- 
definite.  In  order  to  safeguard  against  this  error,  one 
should  perform  a  first-order  analysis.  If  the  error 
message  "WARNING  !  ! !  DECOMPOSITION  FAILED"  shows,  then 
the  stiffness  matrix  is  "singular."  The  input  file 
should  be  checked  for  incorrect  data  of  joint  coordinates 
and  member  connectivity,  which  are  the  common  sources  of 
error . 
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5      VERIFICATION   OF   COMPUTER 
PROGRAM 


5.1    GENERAL 

In  order  to  verify  the  program  [  E  S  C  A  P  ]  ,  the 
following  Sections  will  compare  closed-form  or  published 
solutions  to  the  computer  solutions .  The  second-order 
moments  and  deflections  are  calculated  first  using  the 
closed-form  solution  of  the  differential  equation.  Then 
a  computer  solution  is  generated  and  results  are 
compared . 

5.2   SAMPLE  PROBLEMS   OF  BEAM-COLUMNS 

EXAMPLE  5.2-1   Simply  Supported  Beam-Column  with 

Uniform  Distributed  Load. 
The  simply  supported  beam-column  in  Fig.  15  has  a 
uniform  load  w   of  1  kip/ft  acting  downward.   The  span  L 
=  500  in.  and   I  =  1000  in  4.   Calculate  the  maximum 
deflection  and  moment  at  midspan  based  on  second-order 
elastic  analysis.   Determine  the  elastic  buckling  load 
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and  the  buckled  shape  of  the  member .  Assume  shear  strain 
is  negligible  and  neglect  the  weight  of  the  beam. 

The  maximum  deflection  and  moment  occur  at  midspan. 
The  exact  closed-form  solutions  are  given  by  Timoshenko 
and  Gere1  *  },  and  Chen  and  Lui  [  5  ].  The  general  solutions 
are  given  by  Eqs .  (75)  &  (76). 


y  l12A2sScu_^J__ui]       Eg.  (75) 

5  u4 


where  u^   ;  k  «  Ji  ;  y .-  5wL* 


2        \  EI  °        384EI 


AU  =  MQ    [  2    (seC"  "  1}  ]      Eg.     (76) 


u2 

where  MQ 

= 

wL2 
8 
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(a)     Loading     Diagram 
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(b)        Deflection  Diagram 
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(c)      Moment  Diagram 


Fig.  15     Simply  Supported 
Beam-Column  with  uniform  Load 
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The  computer  results  are  compared  with  these  solutions 
in  TABLES  II  and  III .  These  tables  show  the  computer 
solution  agrees  very  closely  to  the  exact  solution.  The 
amplification  factors  (  AF  )  reported  in  these  tables  are 
the  ratio  of  the  second-order  to  the  first-order  moments 
or  deflections .  Note  that  second-order  effects  are  quite 
significant  at  high  values  of  axial  compressive  force. 

The  theoretical  elastic  critical  load  for  the  beam- 
column  in  Fig.  15  is  given  by  ; 

p    =  %2ei   =  7i2  *  30000  *  1000 
e  '        L2  5002 

=  1184.35  kips.  Eg.    (77) 

The  symbol  P  „  is  used  to  represent  the  elastic  critical 
load. 

The  computer  solution  reports  a  load  factor  LF  = 
2.37,  with  an  applied  axial  force  of  500  kips.  The 
computed  elastic  critical  load  is  therefore  500  *  2 . 37 
=  1185  kips,  versus  the  exact  1184.35  kips.  The  computer 
program  used  a  load  factor  increment  of  . 005  . 

The  computed  relative  joint  displacements  for  the 
buckled  member  are  listed  in  TABLE  IV.  The  buckled  shape 
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is  plotted  in  Fig.  15b.  The  buckled  shape  from  the 
computer  solution  agrees  with  that  predicted  by  the  exact 
solution. 

The  values  of  P  /  P  „  versus  the  moment  and 
deflection  amplification  factors  (MAF,DAF)  are  plotted  in 
Fig.  16.  From  this  figure,  the  following  observations 
are  made  : 

-  the  moment  and  deflection  amplification  factors 
are  non-linear  with  respect  to  P, 

-  the  amplification  factors  are  very  small  for  small 
values  of  P  /  P.  ,on  the  order  of  .05  to  .10, 

-  AF  increases  without  bound  as  P  approaches  P  .. 
These  observations  are  typical  for  second-order  elastic 
response . 


69 


Table  II-  MAXIMUM  BENDING  MOMENTS  FOR 
A  SIMPLY  SUPPORTED  BEAM-COLUMN 

Maximum  Bending  Moments 
(  at  mi ci span  ) 


Axial 
Load 
(  kips  ) 

Exact 
solution 
(  kip-in  ) 

Computer 
solution 
(  kip-in  ) 

MAP 

100 

2851.0 

2851.1 

1.09 

200 

3148.3 

3148.5 

1.21 

300 

3513.1 

3513.2 

1.35 

400 

3971.0 

3971.2 

1.53 

500 

4563.0 

4563.2 

1.75 

Table  III  -  MAXIMUM  DEFLECTIONS  FOR  A  70 

SIMPLY  SUPPORTED  BEAM-COLUMN 


Maximum  Deflections 
(  at  midspan  ) 


Axial 
Load 
(kips  ) 

Exact 
solution 
(  in.  ) 

Computer 
solution 
(  in.  ) 

DAF 

100 

-  2.46965 

-  2.4698 

1.09 

200 

-  2.72135 

-  2.7215 

1.20 

300 

-  3.02999 

-  3.0301 

1.34 

400 

-  3.41735 

-  3.4175 

1.51 

500 

-  3.91793 

-  3.9181 

1.73 

Table  -  IV 


RELATIVE  DISPLACEMENTS  AT  INSTABILITY 
OF  SIMPLY  SUPPORTED  BEAM-COLUMN 
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EXAMPLE  5.2-2   FIXED-FIXED  BEAM-COLUMN  WITH  UNIFORM 

DISTRIBUTED  LOAD 
The  beam-column  shown  in  Fig.  17  is  similar  to  Example 
5.2-1  except  the  end  supports  are  fixed.  We  will 
determine  the  maximum  moment,  deflection,  elastic 
critical  load  and  the  buckled  shape  by  following  the  same 
steps  described  in  Example  5.2-1. 

The  exact  solution  to  the  differential  equation  of 
this  particular  beam-column  is  found  in  Timoshenko  and 
Gere  [  4  },  and  Chen  and  Lui  l  5  ].  The  general  solutions 
are  given  by  Eqs.  (78)  and  (79): 


u2sinu 


where     yJBliX  occurs  at  x  =  —  ;y0  = 


wL4 


2         °        384  EI 


.  _j£  [  3  (tanu  -  u)  ,       BB<  (79) 

12     u2tanu 


Using  w  =  1  kip/ft,  L  =  500  in,  and  I  =  1000  in  4  as 
in  Example  5.2-1,  the  deflection  at  midspan  and  moment  at 
the  end  supports  are  calculated  and  compared  to  the 
computer  solutions  in  TABLES  V  and  VI.  Again,  the 
results  are  in  agreement. 
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(a)  Loading  Diagram 
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(b)  Deflection  Diagram 
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(c)  Moment  Diagram 
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Fig.  17     Fixed-Fixed  Beam-Column 
with  Uniform  Load 


74 
The  theoretical  elastic  critical  load  for  a  fixed- 
fixed  beam-column  is  given  by  : 


*.-**¥■  BQ.    (80) 


m    l^jrj*  30000  .  1000  =  ^  ^ 
5002 


The  computed  load  factor  LF  =  9.4752  with  an  applied 
axial  load  of  500  kips .  The  computed  elastic  critical 
load  is  9.4752  *  500  =  4737.6  kips.  The  load  factor 
increment  was  . 005 .  The  elastic  critical  load  for  a 
fixed-fixed  beam-column  is  four  (4x)  times  larger  than 
the  simply-supported  beam. 

The  computed  relative  joint  displacements  for  the 
buckled  member  are  listed  in  TABLE  VII.  The  buckled 
shape  is  plotted  in  Fig.  17b.  The  buckled  shape  from  the 
computer  solution  agrees  with  that  predicted  by  the  exact 
solution . 

The  relationships  between  P  /  P  „  and  the 
amplification  factors  are  similar  to  those  in  Example 
5.2-1. 
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Table  -  V   MAXIMUM  BENDING  MOMENTS 
FIXED-FIXED  BEAM-COLUMN 


Maximum  Bending  Moments 
(  at  end  supports  ) 


Axial 

Load 

(  kips  ) 

Exact 
Solution 
(  kip-in  ) 

Computer 

Solution 

(  kip-in  ) 

MAF 

100 

1760.6 

1760.7 

1.01 

200 

1786.3 

1786.3 

1.03 

300 

1813.0 

1813.0 

1.04 

400 

1841.0 

1841.0 

1.06 

500 

1870.0 

1870.0 

1.08 

Table  -  VI    MAXIMUM  DEFLECTIONS  76 

FIXED-FIXED  BEAM-COLUMN 


Maximum  Deflections 
(  at  midspan  ) 

Axial 

Load 

(  kips  ) 

Exact 
Solution 
(  in.  ) 

Computer 

Solution 

(  in.  ) 

DAF 

100 

-  .46172 

-  .46174 

1.02 

200 

-  .47176 

-  .47178 

1.04 

300 

-  .48226 

-  .48228 

1.07 

400 

-  .49324 

-  .49326 

1.09 

500 

-  .50474 

-  .50476 

1.12 

Table  -  VII    RELATIVE  DISPLACEMENTS  AT  INSTABILITY 
OF  FIXED-FIXED  BEAM-COLUMN 


RELATIVE  DISPLACEMENTS 
JOINT       D(X)      D(Y)     ROTATION 

1  0  0         0 

2  0-1         0 

3  0         0         0 
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EXAMPLE   5.2-3   HINGED-FIXED  BEAM-COLUMN  WITH  UNIFORM 

DISTRIBUTED  LOAD 
This  problem  is  similar  to  the  previous  two  examples, 
except  this  beam-column  is  restrained  against  rotation  at 
joint  B   and  free  to  rotate  at  joints  A  and  C   (  Fig.  18  ) 
The  maximum  moment,  deflection  at  the  location  of  the 
maximum  moment,  elastic  critical  load,  and  the  buckled 
shape  of  the  member  will  be  determined  by  following  the 
same  procedures  described  in  the  previous  examples . 

The  exact  solution  is  determined  by  superposition. 
The  general  solution  is  a  combination  of  a  simply- 
supported  beam-column  with  a  uniform  load  (  Example  5.2-1 
and  Fig.  19a  )  and  one  with  an  end  moment  acting  at  joint 
B  without  transverse  load  (  Fig.  19b  )  .  For  the  latter 
case  the  general  solutions  are  given  by  Timoshenko  and 
Gere  l  4  ]  : 


=  ^(sin^  _  X)  Eg.    (81) 


y    P   '  sinJcL    L 


Af  =  -Af.  sin^*  Eg.    (82) 

x  D  sinkL 
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Example  5.2-3 


(a)  Load  Diagram 


(b)  Deflection  Diagram 


(c)  Moment  Diagram 


max 


Fig.  18     Hinged-Fixed  Beam-Column 
with  Uniform  Load 
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Example  5.2-3 


tit'fuit 


(a)  uniform  load 


(b)  end  moment 


Fig.  19     Superposition   of 
Lateral  Load  and  Moment 
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The  rotation  at  joint  B  for  the  case  in  Fig.  19b  is  : 


0„  =  ^L  _2_  <_L  -__!_,      EQ.     (83) 
b       3EI  2u     2u        tanu       y 


For  the  casein  Fig.  19a,  the  rotation  at  joint  B  is 


e^^3(tanuru)      *•  «4) 


By  imposing  the  rotational  boundary  condition  at  B,  such 
that  6  b  +   6^  =  0,  the  moment  M  b  is  obtained. 
The  moments  and  the  deflections  at  the  location  of  the 
maximum  moments  for  members  AC  and  CB  are  calculated  and 
compared  to  the  computer  solution  in  TABLES  VI II  thru  X. 
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The  solutions  are  in  agreement. 

The  theoretical  elastic  critical  load  is  l  4  ]  : 


p    -    20.191FJ  _  20.191*30000*1000 
(L)2  (500)2 


Eg.    (85) 
Pe   =  2423  kips 


The  computed  load  factor  LF  =  4 .  8501  for  an  applied  load 
of  500  kips .  The  computed  elastic  critical  load  is 
4.8501  x  500  =  2425.05  kips.  The  load  factor  increment 
was  . 005 .  The  resulting  load  increment  was  . 005  x  500 
kips  =2.5  kips.  The  computed  elastic  critical  load  is 
within  this  tolerance.  Closer  agreement  is  possible  by 
using  a  smaller  load  factor  increment. 

The  computed  relative  joint  displacements  for  the 
buckled  member  are  listed  in  TABLE  XI .  The  buckled  shape 
is  plotted  in  Fig.  18b.  The  buckled  shape  from  the 
computer  solution  agrees  with  that  predicted  by  the  exact 
solution . 
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Table  -  VIII   MAXIMUM  BENDING  MOMENTS 
HINGED-FIXED  BEAM-COLUMN 


Maximum  Bending  Moments 
(  at  Joint  B  of  member  CB  ) 


Axial 

Load 

(  kips  ) 

Exact 
Solution 
(  kip-in  ) 

Computer 

Solution 

(  kip-in  ) 

MAF 

100 

-  2679.4 

-  2679.4 

1.03 

200 

-  2761.0 

-  2761.0 

1.06 

300 

-  2849.9 

-  2849.9 

1.09 

400 

-  2947.0 

-  2947.0 

1.13 

500 

-  3053.7 

-  3053.7 

1.17 

The  deflections  at  the  support  are  0. 
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Table  -  IX   MAXIMUM  MOMENTS  FOR  MEMBER   AC 
HINGED  -  FIXED  BEAM-COLUMN 


Maximum  Bending  Moments 
(  for  member  AC  at  distance  x  from  left  support  ) 


Axial   Distance 
Load     X 
(  kips  )  (in.  ) 


Exact 
Solution 
(  kip-in  ) 


Computer 
Solution 
(  kip-in  ) 


100 

187.11 

1533.2 

1533.3 

200 

186.72 

1608.0 

1608.1 

300 

186.32 

1690.0 

1690.0 

400 

185.91 

1780.2 

1780.3 

500 

185.49 

1880.0 

1880.1 
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Table  -  X   DEFLECTIONS  AT  MEMBER   AC 
HINGED-FIXED  BEAM-COLUMN 


Deflections 
(  for  member  AC  at  distance  x  from  left  support  ) 


Axial   Distance 
Load     X 
(  kips  )  (in.  ) 


Exact 
Solution 
(  in  ) 


Computer 

Solution 

(  in  ) 


DAF 


100 

187.11 

-  .96546 

-  .96636 

1.04 

200 

186.72 

-  1.0083 

-  1.0092 

1.09 

300 

186.32 

-  1.0552 

-  1.0561 

1.14 

400 

185.91 

-  1.1073 

-  1.1076 

1.19 

500 

185.49 

-  1.1646 

-  1.1644 

1.25 
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Table  -  XI    RELATIVE  DISPLACEMENTS  AT  INSTABILITY 

HINGED -FIXED  BEAM- COLUMN 

RELATIVE  DISPLACEMENTS 
JOINT         D(X)        D(Y)         ROTATION 

1  0  0  -  1 

2  0         -  1  .33428 

3  0  0  0 
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EXAMPLE  5.2-4  SECOND-ORDER  ANALYSIS  OF  A  FRAME 
In  previous  examples,  the  given  structures  consisted  of 
a  single-element.   In  this  example,  a  frame  composed  of 
several  elements  will  be  analyzed. 

Since  the  closed-form  solutions  for  frames  are 
difficult  to  derive,  the  published  solution  for  the  given 
frame  is  used  to  verify  the  accuracy  of  the  computer 
solution.  The  effect  of  axial  forces  on  the  behavior  of 
a  rigid  frame  are  illustrated  by  means  moment  diagrams . 

The  rigid  jointed  plane  frame  is  loaded  as  shown  in 
Fig.  20  (  Beaufait  et  al.  )(14]  EI  is  constant  for  all 
members .  The  deflections  and  moments  at  the  joints 
considering  the  effects  of  axial  forces  are  required.  A 
summary  of  the  published  and  computer  solutions  is  shown 
in  TABLES  XII  &  XIII.  The  tables  show  that  the  solutions 
are  in  good  agreement .  The  first  and  second-order 
moments  are  plotted  in  Fig.  21.  From  this  figure,  the 
following  observations  are  made  : 

-  the  largest  amplification  factors  occur  at  the 
support  joints  1  and  4  ,  and  are  on  the  order  of  23%  to 
25%  . 

-  for  each  member  of  the  frame,  the  maximum  moment 
is  located  at  the  end  joint  other  than  the  support 
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joints .  The  amplification  factors  at  the  points  of 
maximum  moments  are  less  than  10%  . 

In  both  examples,  Example  5.2-3  and  this  example, 
the  maximum  moment  do  occur  at  the  member  ends . 


88 


Example  5.2-4 


200  kips 


200  kips 


w  =  -2k/ft 
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^77 
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6  kips 


I5     " 


24  ft 


I  =  250  in 


E  =  30,000  ksi 


Fig.  20        Rigid  Jointed  Frame 
Geometry     and  Loading  Condition 
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Example  5.2-4 


857 


386 

Legends: 

1st  order      (kip-in) 


2nd  order     (kip-in) 


Fig.  21        Moment  Diagrams 
First  and  Second  Order 
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Table  -  XII    SECOND-ORDER  MOMENTS 
OF   A  PLANE  FRAME 


Second-Order  Moments 
Plane  Frame 

Joint 

Published 
Solution 
(  kip-in  ) 

Computer 
Solution 
(  kip-in  ) 

1 

480 

482 

2 

110 

110 

3 

-  610 

-  609 

4 

1153 

1152 

5 

425 

425 

6 

653 

652 
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Table  -  XIII 


SECOND-ORDER  DEFLECTIONS 
OF   A  PLANE  FRAME 


Second-Order  Deflections 
Plane  Frame 


Joint       Published 
Solution 
D  (X)      Rotation 


Computer 
Solution 
D (X)      Rotation 


1 

0 

0 

0 

0 

3 

.8483 

-  .007 

.8526 

-  .007 

4 

.8473 

.0035 

.8495 

.0035 

6 

0 

0 

0 

0 
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EXAMPLE  5.2-5  ELASTIC  BUCKLING  LOAD  OF  A  FRAME 

In  this  example,  the  elastic  critical  load  of  a 
frame  will  be  computed  using  the  method  described  in 
Chapter  4 .  The  results  are  compared  to  the  published 
solution  taken  from  Ghali  and  Neville  [  9  ] .  The  frame 
shown  in  Fig.  22  has  a  load  Q  and  Q/2  applied  at  joints 
2  and  3,  respectively.  Member  2  has  a  flexural  rigidity 
equal  to  EI  ,  while  members  1  and  3  have  a  flexural 
rigidity  equal  to  . 75EI  .  Assuming  that  the  loads  are 
applied  simultaneously,  calculate  the  elastic  critical 
load  Q.  which  causes  buckling  of  the  frame.  Assume  also 
that  the  members  are  made  of  structural  steel  with  the 
following  cross-sectional  and  material  properties: 

I x  =  533  in.  4;  A  ■  19.10  in.  2  ;  E  =  29,000  ksi. 

The  elastic  critical  load  is  determined  by  implementing 
the  procedures  described  in  Section  4.3.  Using  an 
applied  load  Q  =  100  kips,  the  computer  solution  gives  a 
load  factor  LF  =  3 . 94 .  The  computed  elastic  critical 
load  is  therefore  3.94  x  100  =  394  kips.  The  published 
solution  is  .254  EI  /  b  2  =  392.6  kips.  The  error 
between  the  two  solutions  is  very  small  (.36  %  ). 
The  computed  elastic  critical  load  approaches  the 
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exact  solution  when  the  load  factor  increment  is  reduced. 
The  buckled  shape  from  the  computer  solution  agrees  with 
that  predicted  by  the  exact  solution. 

For  the  five  examples  presented,  the  computer 
solutions  are  in  excellent  agreement  with  the  exact 
solutions . 
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Example  5.2-5 
Q  Q/2 
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6b  =  50  ft. 


Qcr  =  .254  El/b 


A  =  19.10  in 


I  =  533  in 


Q  =  100  kips 


Fig.  22  Portal  Frame 

Geometry     and  Loading  Condition 
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EXAMPLE  5.2-6  ELASTIC  SECOND-ORDER  ANALYSIS  OF  A  THREE- 
STORY  FRAME 


In  this  example,  the  capability  of  the  computer  program 
[  E  S  C  A  P  ]  analyzing  larger  structures  is 
demonstrated.  The  elastic  second-order  moments, 
deflections,  and  the  amplification  factors  of  a  three- 
story  structural  steel  frame  will  be  determined  and 
compared  to  the  published  solution. 

The  three-story  frame  in  Fig.  23  l  15 ],  is  subjected 
to  a  combined  gravity  and  lateral  loading  as  shown.  The 
computed  elastic  second-order  moments  and  deflections  are 
shown  in  TABLES  IV  and  XV.  TABLE  IV  shows  that  the  two 
solutions  are  in  agreement. 

The  maximum  lateral  displacement  located  at  the  top 
of  the  structure  (  joint  7  ) ,  is  amplified  about  a  factor 
of  1 .  03  .  This  lateral  displacement  becomes  a  very 
important  design  consideration  in  cases  of  taller 
buildings .  The  lateral  displacement  at  joint  7 
represents  a  story  drift  of  .003H,  where  H  is  the  story 
height  of  the  building.  The  acceptable  lateral 
displacements  for  tall  buildings  are  in  the  range  of 
0.002  to  0.004  of  the  height  of  the  building  '  16  ] . 


Example  5.2-6 
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Fig.  23  Three-Story  Frame 


Geometry  and  Loading 
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Table-XIV  MAXIMUM  SECOND-ORDER  MOMENTS  -  EXAMPLE  5.2-6 


MEMBER       Published 
Solution 
(  kip-in  ) 


1  

2         

3         

4         

5         

6         

7  845 

8  695 

9  579 


Computer 
Solution 
(  kip-in  ) 

MAF 

442 

1.04 

-  26 

0.72 

333 

1.00 

1518 

1.02 

1199 

1.01 

579 

1.01 

846 

1.02 

695 

1.01 

579 

1.01 
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Table-XV  FIRST  and  SECOND-ORDER  LATERAL  DISPLACEMENTS 
EXAMPLE  5.2-6 


JOINT 

First-Order 

2nd-Order 

DAF 

Solution 

Solution 

(  in  ) 

(  in  ) 

1 

0 

0 

0 

2 

0 

0 

0 

3 

0.5419 

0.5603 

1.03 

4 

0.5406 

0.5590 

1.03 

5 

1.1093 

1.1470 

1.03 

6 

1.1080 

1.1460 

1.03 

7 

1.3755 

1.4190 

1.03 

8 

1.3674 

1.4100 

1.03 

CONCLUSIONS 


In  the  analysis  of  structures,  neglecting  the  second- 
order  effects  may  overestimate  the  strength  and  stiffness 
of  a  member  or  frame.  The  elastic  forces  generated 
within  a  member  or  frame  can  be  more  accurately  predicted 
with  the  use  of  an  elastic  second-order  analysis . 

The  second-order  effects  (  amplified  moments  and 
deflections  )  are  of  increasing  importance  as  lighter, 
and  more  flexible  structures  are  constructed.  The  use  of 
higher  strength  materials  and  less  rigid  non-structural 
elements  is  producing  more  flexible  structures  where 
second-order  effects  are  of  greater  importance. 

The  behavior  of  multistory  frames  subjected  to 
either  gravity  loads  only  or  under  combined  gravity  and 
lateral  loading  cannot  be  accurately  predicted  by  a 
first-order  elastic  analysis,  when  the  compressive  axial 
load  level  is  significant  '  17  ] .  The  lateral 
displacements  may  be  rather  large  and  may  not  be  within 
the  acceptable  range  of  lateral  displacements  for  tall 
buildings .  An  elastic  second-order  analysis  such  as 
[  E  S  C  A  P  ]  can  be  used  to  provide  an  accurate 
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100 
estimate  of  the  lateral  displacements  and  the  elastic 
internal  forces  generated  in  the  structure.  If  the 
lateral  displacements  become  too  large,  they  may  become 
noticeable  and  disturbing  to  the  building  occupants  [  16 

When  second-order  effects  are  included  in  the 
analysis,  a  more  rational  design  of  individual  members 
may  be  possible.  The  beam-column  interaction  equation  in 
LRFD  is  based  on  the  use  of  elastic  second-order  moments . 
LRFD  permits  direct  application  of  computer  generated 
elastic  second-order  solutions,  such  as  those  generated 
by  the  computer  program  [  E  S  C  A  P  ] . 

Structures  with  very  slender  members  may  fail  by 
elastic  instability.  For  such  structures,  elastic 
stability  analysis  is  directly  applicable  when 
determining  the  elastic  buckling  load. 

Typical  members  or  frames  in  structural  steel 
buildings  will  fail  by  inelastic  instability,  rather  than 
elastic  instability.  In  such  cases,  elastic  frame 
stability  analysis  is  not  directly  applicable  for 
predicting  the  structural  strength.  However,  the  LRFD 
beam-column  interaction  equations,  which  predict  frame 
instability,  require  the  use  of  effective  length  factors. 
An  elastic  frame  stability  analysis  provides  a  powerful 
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tool  to  determine  the  effective  length  factors. 

The  program  [  E  S  C  A  P  ]  was  written  to  provide 
exact  elastic  second-order  analysis  based  on  the  use  of 
stability  functions .  The  program  was  verified  by 
comparison  with  published  closed-form  solutions . 

The  program  [  E  S  C  A  P  ]  also  permits  elastic 
stability  analysis  of  members  and  plane  frames . 
Instability  is  determined  by  increasing  the  applied  load 
increments  until  the  elastic  second-order  stiffness 
matrix  becomes  singular.  Rather  than  checking  the 
singularity  by  computing  the  determinant,  singularity  was 
detected  by  checking  for  a  zero  or  negative  diagonal 
element  using  Cholesky's  decomposition  of  the  stiffness 
matrix.  This  provides  a  simple  technique  to  incorporate 
stability  analysis  into  an  elastic  second-order  program. 

The  accuracy  of  the  program  [ESCAP]  in 
computing  elastic  stability  loads  was  verified  by 
comparison  with  known  closed-form  solutions .  The 
solutions  are  in  excellent  agreement. 

The  future  of  structural  design  practice  will 
include  second-order  effects  in  the  computer  analysis  and 
design  of  structures .  Structural  analysis  programs  like 
[ESCAP  ]  will  become  routine. 


APPENDIX 
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APPENDIX   A 

USERS  MANUAL 
ELASTIC  SECOND-ORDER  COMPUTER  ANALYSIS 
PROGRAM 

Al   INTRODUCTION 

The  Elastic  Second-Order  Computer  Analysis  Program  [  E 
S  C  A  P  ]  was  developed  to  provide  designers  a  tool  in  the 
analysis  and  design  of  beam-columns  and  frames .  The 
program  [  E  S  C  A  P  ]  has  the  capability  of  performing 
first-order  linear  elastic  analysis  as  well  as  elastic 
second-order  analysis .  The  program  was  also  designed  to 
perform  elastic  stability  analysis  of  frames  and  beam- 
columns  .  The  program  can  be  used  to  perform  structural 
analysis  on  any  two-dimensional  structures  subjected  to  in- 
plane  bending. 

The  program  can  be  loaded  in  any  IBM  PC  or  compatible 
with  a  hard-disk  or  floppy.  It  is  written  in  single- 
precision  FORTRAN  with  the  IBM  Professional  Fortran 
Compiler. 
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A2   CAPABILITY 

The  program  E  S  C  A  P  has  the  following  capabilities  : 

a.  Performs  first-order  linear  analysis  of 
plane  frames,  beams,  columns,  plane  trusses 
and  beam-columns . 

b.  Performs  elastic  second-order  analysis. 

c.  Performs  elastic  stability  analysis. 

d.  Determines  the  location  and  magnitude  of  the 
maximum  elastic  second-order  moment  within 
each  member. 

e.  Has  the  capability  of  analyzing  multiple  load 
patterns . 


A3   PROGRAM  PARAMETERS 

The  following  parameters  were  used  in  the  computer  program 
[  E  S  C  A  P  ]  : 

PID  -  Problem  Identification.  This  is  the  name  of  the 
input  file  with  a  brief  description  of  the  problem.  It 
consists  of  12  characters  per  line  with  a  maximum  of  three 
(3)  lines  provided  in  the  input  file. 

TITLE  is  a  heading  that  describes  the  input 
parameters  in  a  read  statement. 
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NJ  -  number  of  joints;  Maximum  of  40.   The  user  has  an 

option  to  increase  the  size  by  increasing  the  dimension  of 

the  joint  coordinates  X  (  )  and  Y  (  ) . 

NDELTA  -  is  a  code  that  determines  the  type  of 

analysis . 

0  :  1st  order  analysis 

1  :  2nd  order  analysis 

NCYCLE  -  it  is  the  maximum  number  of  iterations.  The 
program  will  terminate  if  after  NCYCLE  the  second-order 
solution  has  not  converged,  and  will  print  the  output  of 
the  last  computed  values  of  member  forces  and  joint 
deflections.  It  normally  takes  between  4  to  10  iterations 
before  the  solution  converges . 

EPS1  -  is  the  tolerance  factor  used  to  test 
convergence . (  recommend   .005  ) 

STEP  -  is  the  load  factor  increment  used  in  the 
elastic  stability  analysis. (  recommend   .005  ) 

LB  -  is  the  code  that  determines  if  buckling  load 
analysis  is  desired. 

0  :  NO 

1  :  YES 

X  (  J  )  -  is  the  x  -  coordinate  of  a  joint  J,  in  feet. 
Y  (  K  )  -  is  the  y  -  coordinate  of  a  joint  K  in  feet. 
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MCL  -  is  the  joint  release  code. 

0  :  joint  is  free 

1  :  joint  is  fixed. 

NDOF  -  is  the  total  number  of  degrees  of  freedom. 
(  equal  to  3  *  NJ  ) 

NRJ  -  number  of  restrained  joints. 

S  (  I,  J  )  -  is  the  stiffness  matrix  of  the  structure 
or  a  member 

NRL  -  number  of  restrained  degrees  of  freedom  (DOF's) 
that  are  fixed  and  have  0  displacements. 

MNDOF  =  NDOF  -  NRL  ,  is  the  total  number  of  degrees  of 
freedom  that  are  free  to  rotate  and  translate. 

NM  -  is  the  number  of  members.   Maximum  is  80. 
The  user  has  the  option  to  increase  the  size  by  changing 
the  dimension  of  the  other  parameters  in  the  main  program. 

MPC  -  Matrix  Print  Code  for  the  stiffness  matrix. 

0  :  suppress  output 

1  :  activate  output 

GX  -  is  the  unit  weight  of  the  material;  490  pcf  for 
structural  steel . 

IHINGE  (  I  )  -  is  a  joint  code  that  describes  the 
type  of  restraints  imposed  at  the  joints .   It  is  used  to 
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modify  the  stiffness  of  a  member  that  is  pinned  at  one  end. 
For  members  without  supports  : 

0  :  pinned  at  start  joint  and  rigid  at  the 

end  joint 

1  :  rigid  joints  at  both  ends 

2  :  joints  pinned  at  both  ends 

For  members  with  support  joints  (hinged  or  fixed) : 
1  :  if  the  beginning  joint  is  a  support 
joint 

0  :  if  the  beginning  joint  is  pinned  and  not 

a  support  joint 
JFIX  (  J  )  ,  KFIX  (  K  )  -   are  joint  codes  used  to 
modify  the  load  vector  and  the  force  vector  for  a  pin 
joint . 

1  :  for  a  fixed  support  or  rigid  joint 

0  :  for  a  hinged  support  or  pinned  joint 
JM  (  I  )  -  beginning  joint  of  a  member 
KM  (  I  )  -  end  joint  of  a  member 
AM  (  I  )  -  area  of  cross-section,  sq.  in. 
ZM  (  I  )  -  moment  of  inertia;  in  4 
EM  (  I  )  -  modulus  of  elasticity  of  the  material 

Example  :  Enter  29  (ksi)  for  steel 
SPAN  -  member  length  in  inches 
SM  (  I  )  -  member  length  in  feet 
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CX  (  I  )  -  cosine 
CY  (  I  )  -  sine 

WT  -  total  weight  of  structure  in  tons. 
INODE  (  )  -  is  a  degree  of  freedom  (  DOF  )  at  the  ends 
of  a  member  : 

(  1  )  -  DOF  at  joint  JM  (  I  )  along  the  x  -  axis 

(  2  )  -  DOF  at  joint  JM  (  I  )  along  the  y  -  axis 
(  3  )  -  rotation  at  joint  JM  (  I  ) 
(  4  )  -  DOF  at  joint  KM  (  I  )  along  the  x  -  axis 
(  5  )  -  DOF  at  joint  KM  (  I  )  along  the  y  -  axis 
(  6  )  -  rotation  at  joint  KM  (  I  ) 

Q    (  I, J  )  -  stiffness  matrix  for  the  restrained 

DOF's 
NLP  -  number  of  loading  patterns 
LP  -  load  pattern 
NLJ  -  number  of  loaded  joints 
WX  -  joint  load  in  x  -  direction  (  kips  ) 
WY  -  joint  load  in  y  -  direction  (  kips  ) 
WM  -  joint  moment  (  kip- ft  ) 

W  (  )  -  sum  of  joint  loads  in  x  ,  y  or  moments 
FEA  (  )  -  fixed-end  actions  or  forces  in  members  due 
to  externally  applied  loads . 

MEA  (  )  -  member  end  actions  after  the  analysis 
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NMC  -  number  of  members  with  concentrated  loads 
NCL  -  number  of  concentrated  loads  on  a  member  ( 
limited  to  4) .  The  user  can  increase  the  size  by  changing 
the  dimension  of  LCN  in  the  main  program.  The  location  of 
loading  is  expressed  as  a  fraction  of  the  span  measured 
from  the  left  support. 

FX  -  concentrated  load  in  x  -  direction  (  kips  ) 
FY  -  concentrated  load  in  y  -  direction  (  kips  ) 
MO  -  moment  in  z  -  direction  (  kip-ft  ) 
NMD  -  number  of  members  with  uniform  loads 
NUL  -  number  of  uniform  loads  on  a  member 
SUMWY  (  )  -  sum  of  transformed  uniform  distributed 
loads  on  the  member  along  the  y  -  axis  . 
MWC  -  member  weight  code. 

0  :  exclude  members  weight  in  the  analysis 

1  :  include  member  weight  in  the  analysis 
R  (  )  -  reaction  vector 

6  (  )  -  modified  fixed  end  actions 
RAL  (  )  -  identifies  the  restrained  DOF 
D  (   )  -  deflection  vector 
XNEW  (  )  -  new  x-coordinates 
YNEW  (  )  -  new  y-  coordinates 

Dl,  D2,  D3,  D4,  D5,  D6  -  deflections  at  member's  local 
axis 
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FIX  (  )  -  calculated  fixed-end  actions  during  second- 
order  iteration  analysis 

ST  (  )  -  modified  stiffness  matrix  for  second-order 
effects . 

LF  -  load  factor  use  in  load  buckling  analysis.  It  is 
equal  to  1  at  the  start  of  the  analysis. 

ITER  -  counts  number  of  times  the  load  factor 
increment  STEP  was  added  to  LF  before  the  buckling  load  is 
achieved . 

QT  -  restrained  stiffness  matrix  for  the  second-order 
analysis 

P  -  is  the  axial  load  or  force 

XKL  =  sqrt  (  P  /  EI  ) 

VA,  VB  -  vertical  reactions  considering  axial  load 
effects 

MA,  MB  -  moment  reactions  considering  axial  load 
effects 

COMP  -  stability  functions  for  axial  compression 

TENS  -  stability  functions  for  axial  tension 

FLEXMAX  -  maximum  moment  anywhere  in  the  span 

YMAX  -  deflection  of  FLEXMAX.  (Note:  this  may  or  may 
not  be  the  maximum  deflection  ) 

XMAX  -  location  of  FLEXMAX  and  YMAX 
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A4  LIST  OF  SUBROUTINES 

a.  E  S  C  A  P  is  the  main  program.  It  reads  the 
input  file  and  performs  the  first-order  analysis. 

b.  NODE  is  a  subroutine  that  generates  the  degrees  of 
freedom  for  each  node. 

c.  STIFF  calculates  and  assembles  the  global  first- 
order  stiffness  matrix. 

d.  MATPRT  prints  the  global  stiffness  matrix. 

e.  DECOMP  is  the  CHOLESKY's  matrix  decomposition 
program.  It  decomposes  the  global  stiffness  matrix  and 
stores  the  matrix  in  S  (  I, J  )  or  ST  (  I,  J  ) . 

f.  LOAD  modifies  the  fixed-end  actions  FEA  for  a 
member  that  is  pinned  at  the  beginning  node  or  both. 

g.  SOLVER  solves  the  equations  for  unknown 
deflections  and  stores  them  at  vector  D  (  ) . 

h.  MFORCE  calculates  the  first-order  member  end 
actions  in  the  local  axis  and  stores  them  at  vector  MEA 

i .   PNDELTA  performs  the  second-order  analysis . 

j.  XDELTA  generates  the  global  second-order  stiffness 
matrix  and  stores  the  matrix  at  vector  ST  . 

k.  DFORCE  solves  the  members  end  actions  (MEA)  to 
include  second-order  effects . 

1.    CALMAX  determines  the  location  of  the  maximum 
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moment  for  a  uniformly  loaded  member  with  end  moments .  It 
also  calculates  the  elastic  second-order  moments  and 
deflections. 


A5  INPUT  DATA  PREPARATION 

An  editor  (such  as  PE  FORTRAN  )  is  used  to  prepare  the 
input  file.   The  parameters  of  every  read  statements  are 
entered  in  separate  lines  and  are  format  free.   A  sample 
copy  of  the  input  file  is  shown  in  Appendix  B. 
PROCEDURES : 

A.  READ  PID  -  in  alphanumeric  characters  and  consists 
of  12  characters  per  line  with  a  maximum  of  3  lines . 

B.  READ    TITLE    -    Enter    >    number    of 
joints , NDELTA ,NCY CLE,  EPS1,STEP,LB 

C.  READ  NJ,  NDELTA,  NCYCLE,  EPS1 

D.  READ  STEP,  LB 

E.  READ  TITLE  -  Enter>  Joint  Coordinates 

F.  READ  J,  X  (  J  ) ,  Y  (  J  ) 

J  :  is  the  joint  or  node 

G.  READ  TITLE  -  Enter>  Number  of  restrained  joints 
H.  READ  NRJ 

I .  READ  TITLE  -  Enter>  Boundary  Conditions   (Note : 
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these  are  the  boundery  conditions  of  the  restrained  joints 
in  procedure  G  &  H) 

J.  READ  J,  ,MCL  (  1  ),  MCL  (  2  )  ,  MCL  (  3  ) 
J  :  is  the  joint  or  node 
MCL  (  1  )  :  DOF  in   X  -  DIRECTION 
MCL  (  2  )  :  DOF  in   Y  -  DIRECTION 
MCL  (  3  )  :  DOF  in  ROTATION 
K.  READ  TITLE  -  Enter>  Number  of  members , print 

code 
L.  READ  NM,MPC 

M.  READ  TITLE  -  Enter>  Member  Data 
N.  READ  I,J,K,A,Z,E,GX 

O.  READ  IHINGE  (I),JFIX  (J),  KFIX  (K) 
I  :  is  the  member 
J  :  beginning  node 
K  :  ending  node 

A  :  area  of  cross-section  (  in  2  ) 
Z  :  moment  of  inertia  (  in  *  ) 
E  :  section  modulus  (  KSI/1000  ) 
GX  :  unit  weight  of  member  (  LBS/FT  3  ) 
P.  READ  TITLE  -  Enter>  Number  of  Load  Patterns 
Q.  READ  LP 

R.  READ  TITLE  -  Enter>  Number  of  loaded  joints 
S.  READ  NLJ  ;  if  NLJ  is  equal  to   0   go  to 
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step  V  . 
T.  READ  TITLE  -  Enter>  Loaded  joints 
U.  READ  J,WX,WY,WM 

J  :  is  the  joint  or  node 
V.  READ  TITLE  -  Enter>Number  of  members  with 

concentrated  loads 
W.  READ  NMC  ;  if  NMC  is  equal  to   0  ,  go  to 

step  AA 
X.  READ  TITLE  -  Enter>  Members  concentrated  loads 
Y.  READ   M,  NCL 
Z.  READ   K,  A,  FX,  FY,  MO 
K  :  load  orientation 

0  :  member  axis 

1  :  structure  or  global  axis 

A  :  location  of  loading  (  fraction  of  the 
span  ) 
AA.   READ  TITLE  -  Enter>  Number  of  members  with 
distributed  loads 

BB.  READ  NMU  ;  if  NMU  is  equal  to   0   go 

to  step   FF 
CC.  READ  TITLE  -  Enter>  Members  uniform  loads 
DD.  READ   M,NUL 
EE.  READ   K,FX,FY,MO 

M  :  name  of  member 
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K  :  load  orientation 

0  :  local  or  member  axis 

1  :  global  or  structure  axis 

MO  :  is  expressed  in  kip-ft  /  ft 
FF.  READ  TITLE  -  Enter>  Member  weight  code 
GG.  READ   MWC 

0  :  exclude  member  weight 

1  :  include  member  weight 

A6  OUTPUT  DATA 

The  computer  output  shows  the  joint  coordinates,  members 
data,  joint  loading,  and  the  member  loading  in  a  load 
pattern.  It  also  prints  the  joint  deflections,  member  end 
actions  and  the  reaction  forces .  The  location  of  the 
maximum  moment,  its  magnitude  and  deflection  are  also 
shown.  The  load  factor  LF  is  printed  and  used  in 
determining  the  actual  buckling  load,  and  is  obtained  by 
multiplying  the  load  factor  LF  with  the  member' s  load  or 
joint  load.   A  sample  output  data  is  shown  in  Appendix  B. 
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A7   RUNNING  THE  PROGRAM 

The  following  steps  describe  how  to  run  the  Elastic  Second- 
Order  Computer  Analysis  Program  [  E  S  C  A  P  ] : 

1.  Name  the  input  file  described  in  Section  A6  above. 
The  input  file  is  named  when  saving  the  input  data. 

2 .  Go  to  the  Fortran  Directory . 

3 .  To  run  the  program  : 

ENTER  >  XFRAME 
RETURN 

4 .  The  program  will  prompt > 

What  is  the  name  of  the  input  file? 

5.  ENTER  >  Name  of  input  file.   Make  sure  the  input 
file  is  saved  in  the  Fortran  directory. 

6.  RETURN 

7 .  Execution  of  the  program  is  completed  when  the 
message  Execution  Terminated  :  0  shows  in  the  screen. 

8.  The  output  file  is  automatically  sent  to  the  file 
named  SOLOUT  in  the  Fortran  directory. 


APPENDIX      B 


COMPUTER   INPUT  AND   OUTPUT   DATA 
For  Example   5.2-4 
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INPUT  FILE  -  SAMPLE. IN 
RIGID  FRAME 
EXAMPLE  5.2-4 
NUMBER  OF  JOINTS, NDELTA, NCYCLE, EPS1,  STEP,  LB 
6,1,1000, .005 
.005,0 

JOINT  COORDINATES 
1,0,0 
2,0,12 
3,0,20 
4,24,20 
5,24,6 
6,24,0 

NUMBER  OF  RESTRAINED  JOINTS 
2 
BOUNDERY  CONDITIONS 

1,1,1,1 

6,1,1,1 

NM,MPC,MCODE 

5,0,1 

MEMBER  DATA 

1,1,2,13.3,250,30,490 

1,1,1 

2,2,3,13.3,250,30,490 

1,1,1 

3,3,4,13.3,250,30,490 

1,1,1 

4,4,5,13.3,250,30,490 

1,1,1 
5,5,6,13.3,250,30,490 

1,1,1 

LOAD  PATTERNS 

1 

NUMBER  OF  LOADED  JOINTS 

4 

LOADED  JOINTS 

2,12,0,0 

3,0,-200,0 

4,0,-200,0 

5,-6,0,0 

NUMBER  OF  MEMBERS  W/  CONCENTRATED  LOADS 

0 

NUMBER  OF  MEMBERS  WITH  DISTRIBUTED  LOADS 

1 

LOADED  MEMBER 

3,1 

1,0,-2,0 
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MEMBERS  WEIGHT  CODE 
0 
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"  OUTPUT  DATA  " 

ELASTIC   SECOND-ORDER  COMPUTER 
ANALYSIS  OF  BEAM-COLUMNS  AND  FRAMES 
PLANE  FRAME  ANALYSIS  BY  STIFFNESS  METHOD 

***************** 

***   GENE    MUNOZ   *** 
***************** 


****** 

*  J  O  I  N  T 
****** 


***** 

DATA* 
***** 


JOINT      COORDINATES 
X        Y 
(FEET)    (FEET) 


RESTRAINT 
CONDITION 
X   Y   R 


0.000    0.000 


2 

0.000 

12.000 

3 

0.000 

20.000 

4 

24.000 

20.000 

5 

24.000 

6.000 

6 

24.000 

0.000 

******   ****** 
*MEMBER   DATA* 


******* 


*  *  *  *  * 


EMBE 

R    LENGTH 

CROSS-SECTION 

PROPERTIES 

MODULUS  OF 

AREA 

I(ZZ) 

ELASTICITY 

(FEET) 

(INS**2) 

(INS**4) 

(KSI/1000) 

1 

12.000 

13.300 

250.00 

30.000 

2 

8.000 

13.300 

250.00 

30.000 

3 

24.000 

13.300 

250.00 

30.000 

4 

14.000 

13.300 

250.00 

30.000 

5 

6.000 

13.300 

250.00 

30.000 
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TOTAL  WEIGHT  OF  FRAME: 


1.448    TONS 


******************** 
*  *  *  FIRST  ORDER  ***** 
****  ANALYSIS  ****** 
******************** 

LOADING   PATTERN       1    OF       1    LOADING   PATTERNS 


*************** 

*JOINT    LOADING* 
*************** 


JOINT 


2 
3 
4 
5 


LOADING 

F(X) 

F(Y) 

COUPLE 

(KIPS) 

(KIPS) 

(KIP-FEET) 

12.000 

0.000 

0.000 

0.000 

-200.000 

0.000 

0.000 

-200.000 

0.000 

-6.000 

0.000 

0.000 

********   ******   * 
♦MEMBER    LOADING* 


****** 


*  * 


NO  MEMBERS  WITH  CONCENTRATED  LOADS  IN  PATTERN 


MEMBER  UNIFORMLY 

ORIENTATION   W(X) 
(AXES)     (KIPS/FT) 


DISTRIBUTED   LOADING 

W(Y)  M(O) 

(KIPS/FT)      (K-FT/FT) 


0.000 


-2.000 


0.000 
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********************** 

*  MEMBER  WEIGHTS  NOT  INCLUDED  IN  ANALYSIS  * 
********************** 

LOADING  PATTERN   1  OF   1  LOADING  PATTERNS 


*********         ********* 

*JOINT         DEFLECTIONS* 
******************* 


JOINT  DEFLECTIONS 

D (X)  D (Y)  ROTATION 

(INCHES)  (INCHES)  (RADS) 


1  0.0000E-01      0.0000E-01      0.0000E-01 

2  3.0288E-01    -8.0369E-02    -2.5911E-03 

3  6.4644E-01    -1.3395E-01    -6.3060E-03 

4  6.4031E-01    -1.3553E-01      3.8776E-03 

5  1.6280E-01    -4.0658E-02    -4.2348E-03 

6  0.0000E-01   0.0000E-01   0.0000E-01 


LOADING  PATTERN   1  OF   1  LOADING  PATTERNS 


******************* 

♦MEMBER      END      ACTIONS* 
******************* 


MEMBER 
1 
2 


END 
JOINTS 

1 
2 

2 

3 


X-FORCE 
(KIPS) 

2.2269E+02 
-2.2269E+02 

2.2269E+02 
-2.2269E+02 


Y-FORCE 
(KIPS) 

3.5060E+00 
-3.5060E+00 

-8.4940E+00 
8.4940E+00 


MOMENT 
(K-IN) 

3.8739E+02 
1.1748E+02 

-1.1748E+02 
-6.9794E+02 


8.4941E+00 
8.4941E+00 


2.2688E+01 
2.5312E+01 


6.9794E+02 
-1.0757E+03 


123 


4 
5 


2.2531E+02 
-2.2531E+02 


8.4941E+00 
•8.4941E+00 


1.0757E+03 
3.5134E+02 


5 
6 


2.2531E+02 
■2.2531E+02 


2.4941E+00 
•2.4941E+00 


-3.5134E+02 
5.3091E+02 


LOADING  PATTERN   1  OF   1  LOADING  PATTERNS 


********      ******** 

♦REACTIVE      FORCES* 
*********      ******* 


JOINT 


REACTIONS 
R(X)  R(Y)  MOMENT 

(KIPS)  (KIPS)  (KIP-INCH) 


1 

-3.5060E+00 

2.2269E+02 

3.8739E+02 

6 

-2.4  941E+00 

2.2531E+02 

5.3091E+02 
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***************** 
*  SECOND  ORDER  *** 
**         ANALYSIS  *** 

***************** 


LOADING   PATTERN      1    OF      1    LOADING   PATTERNS 


*********         ********* 

*    J   O   I    N   T         DEFLECTIONS* 
******************* 


JOINT  DEFLECTIONS 

D (X)  D (Y)  ROTATION 

(INCHES)  (INCHES)  (RADS) 

1  0.0000E-01   0.0000E-01   0.0000E-01 

2  4.0855E-01  -8.0162E-02  -3.7614E-03 

3  8.5261E-01  -1.3360E-01  -7.0316E-03 

4  8.4646E-01  -1.3587E-01   3.5363E-03 

5  2.0170E-01  -4.0761E-02  -5.2386E-03 

6  0.0000E-01   0.0000E-01   0.0000E-01 

CONVERGENCE  OCCURED  AT  ITER   = 
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LOADING  PATTERN   1  OF   1  LOADING  PATTERNS 


******************* 

*MEMBER      END      ACTIONS* 
******************* 


MEMBER 

END 

X-FORCE 

Y-FORCE 

MOMENT 

JOINTS 

(KIPS) 

(KIPS) 

(K-IN) 

1 

1 

2.2212E+02 

3.4786E+00 

4.8161E+02 

2 

-2.2212E+02 

-3.4786E+00 

1.1006E+02 

2 

2 

2.2212E+02 

-8.5214E+00 

-1.1006E+02 

3 

-2.2212E+02 

8.5214E+00 

-6.0936E+02 

3 

3 

8.5212E+00 

2.2115E+01 

6.0936E+02 

4 

-8.5212E+00 

2.5885E+01 

-1.1522E+03 

4 

4 

2.2588E+02 

8.5213E+00 

1.1522E+03 

5 

-2.2588E+02 

-8.5213E+00 

4.2502E+02 

5 

5 

2.2588E+02 

2.5213E+00 

-4.2502E+02 

6 

-2.2588E+02 

-2.5213E+00 

6.5212E+02 

LOADING   PATTERN      1    OF      1    LOADING   PATTERNS 


********       ******** 

♦REACTIVE      FORCES* 
*********      ******* 


JOINT 


REACTIONS 
R(X)  R(Y)  MOMENT 

(KIPS)  (KIPS)  (KIP-INCH) 


1 
6 


-3.4786E+00 


-2.5213E+00 


2.2212E+02 


2.2588E+02 


4.8161E+02 


6.5211E+02 
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LOADING  PATTERN   1  OF   1  LOADING  PATTERNS 


***************** 

♦MAXIMUM    MOMENTS* 
***************** 


MEMBER  MOMENTS 

XMAX         YMAX         MOMENT 
(INCH)        (INCH)       (KIP-INCH) 

1  0.0000E-01 

0.0000E-01 

4.8161E+02 

2  9.5948E+01 

-1.3360E-01 

-6.0936E+02 

3  2.8799E+02 

-1.3587E-01 

-1.1522E+03 

4  0.0000E-01 

-1.3587E-01 

1.1522E+03 

5  7.1960E+01 

0.0000E-01 

6.5212E+02 


ANALYSIS  COMPLETED 
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